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Materials  Research  Laboratory  is  the  Program  Manager  and  Dr.  S.  D. 
Manning  is  the  Principal  Investigator.  This  is  Phase  I  of  a  three  phase 
program. 


A  durability  analysis  methodology  has  been  developed  under  Phase  I 
for  satisfying  the  Air  Force's  durability  design  requirements  for  advanced 
metallic  airframes.  This  effort  is  documented  in  five  separate  volumes: 

Volume  I  -  Durability  Methods  Development  -  Phase  I  Summary 

Volume  II  -  Durability  Analysis:  State-of-the-Art  Assessment 

Volume  III  -  Structural  Durability  Survey:  State-of-the-Art 

Assessment 

Volume  IV  -  Initial  Quality  Representation 

Volume  V  -  Durability  Analysis  Methodology  Development 


The  Phase  I  effort  is  summarized  in  Volume  I.  This  report  (Volume 
V)  documents  in  detail  the  durability  analysis  concepts  and  procedures 
developed  under  Phase  I. 
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DEFINITION  OF  TERMS 


1.  Deterministic  Crack  Growth  -  Crack  growth  parameters  are  treated  as 

deterministic  values  resulting  in  a  single  value  prediction  for  crack 
length. 

2.  Durability  -  "The  ability  of  the  airframe  to  resist  cracking  (includ¬ 
ing  stress  corrosion  and  hydrogen  induced  cracking),  corrosion,  thermal 
degradation,  delamination,  wear,  and  the  effects  of  foreign  object 
damage  for  a  specified  period  of  time"  [1] . 

3.  Durability  Analysis  -  An  analysis  for  quantifying  the  extent  of  struc¬ 
tural  damage  as  a  function  of  service  hours  used  for  ensuring  compliance 
with  the  Air  Force’s  durability  design  requirements  for  advanced  metallic 
airframes. 

4.  Economic  Life  -  "...  occurrence  of  widespread  damage  which  is  uneconom¬ 
ical  to  repair  and,  if  not  repaired,  could  cause  functional  problems 
affecting  operational  readiness"  [1] . 

5.  Economic  Life  Criteria  -  Analytical  guidelines  for  quantifying  economic 
life. 

6.  Economic  Repair  Limit  -  Maximum  damage  size  that  can  be  economically 
repaired  (e.g.,  repair  0.030"-. 050"  radial  crack  in  fastener  holes  by 
reaming  hole  to  next  size) . 

7.  Equivalent  Initial  Flaw  Size  (EIFS)  -  A  hypothetical  crack  assumed  to 
exist  in  the  structure  prior  to  service.  It  characterizes  the  equiva¬ 
lent  effect  of  actual  initial  flaws  in  a  structure  detail.  Values  are 
usually  determined  by  back  extrapolating  fractographic  results. 

8.  Initial  Quality  -  "A  measure  of  the  condition  of  the  airframe  relative 
to  flaws,  defects,  or  other  discrepancies  in  the  basic  materials  as 
introduced  during  manufacture  of  the  airframe"  [1] . 

9.  Probability  of  Crack  Exceedance  -  Refers  to  the  probability  of  exceeding 
a  specified  crack  size  at  a  given  service  time.  It  can  be  determined 
from  the  statistical  distribution  of  crack  sizes  and  it  is  the  funda¬ 
mental  quantity  for  predicting  economic  life. 

10.  Time-To-Crack-Initiation  (TTCI)  -  The  time  or  service  hours  required  to 
initiate  a  specified  (observable)  crack  size  in  a  structural  detail. 
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SECTION  I 


INTRODUCTION 


Current  U.  S.  Air  Force  structural  integrity  and  durability  design 
specifications  [l-3]  require  that  airframe  components  be  designed  such 
that  the  economic  life  be  analytically  predicted.  The  conventional 
fatigue  analysis,  while  capable  of  estimating  design  life,  does  not  lend 
itself  to  predicting  the  economic  life,  nor  is  it  capable  of  providing 
a  definition  of  economic  life.  In  predicting  the  economic  life,  the  entire 
crack  population  should  be  taken  into  account,  and  hence  the  statistical 
approach  is  essential.  Literature  describing  the  current  practice  re¬ 
garding  crack  growth  damage  accumulation  with  applications  to  structural 
safety,  durability,  damage  tolerance,  and  reliability  is  available  [4-9]  . 

In  particular,  the  durability  and  the  economic  life  requirements  have  been 
reviewed  in  detail  in  Ref.  9. 

The  cost  of  maintenance,  including  the  costs  of  inspection,  repair, 
rework,  replacement,  down  time,  etc.,  for  aircraft  structures,  in  order 
to  fulfill  the  requirements  of  safety,  durability,  danage  tolerance  and 
reliability,  is  of  practical  importance.  It  is  well-known  that  after  a 
certain  service  life,  referred  to  as  the  economic  life,  either  the  cost 
of  maintenance  or  the  number  of  cracks  exceeding  the  economical  repair 
crack  size  increases  so  rapidly  that  the  durability  requirement  cannot 
be  satisfied.  In  an  attempt  to  demonstrate  statistically  the  existence 
of  the  economic  life  quantitatively,  an  exploratory  investigation  was 
made  to  estimate  the  service  cracks  and  maintenance  cost  [10],  In 
Ref.  10,  however,  restrictive  limitations  were  imposed,  for  simplicity, 
in  the  exploratory  studies. 

In  this  volume,  an  analytical  durability  methodology  that  is  capable 
of  quantitatively  predicting  the  economic  life  of  advanced  metallic  air¬ 
craft  structures  has  been  developed.  The  methodology  is  based  on  sound 
analytical  and  statistical  approaches.  It  accounts  for  various  service 
conditions,  such  as  any  type  of  initial  fatigue  quality,  crack  growth 
damage  accumulation,  loading  spectra,  material/structural  properties, 
usage  change,  inspection  and  repair  maintenance,  etc.  Hence,  the  method- 
ology  presented  herein  is  general  enough  for  practical  applications.  The 
formulation  allows  for  the  determination  of  the  economic  life  using  either 
one  of  the  following  two  criteria;  (i)  a  rapid  increase  of  the  number  of 
crack  damages  exceeding  the  economic  repair  crack  size,  and  (ii)  rapid 
increase  of  the  maintenance  cost.  The  economic  repair  crack  size  ae  is 
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defined  as  the  crack  size  below  which  the  least  expensive  repair  procedure 
can  be  used,  such  as  remaing  the  fastener  holes  to  the  next  hole  size.  The 
size  ag  is  usually  between  0.03"  and  0.05"  depending  on  the  location  and 
the  fastener  hole  size. 

The  durability  critical  component  may  also  be  subjected  to  a  sched¬ 
uled  (nonperiodic)  inspection  and  repair  maintenance  procedure  as  shown 
in  Fig.  (1),  in  which  a  component  without  a  maintenance  procedure  is  a 
special  case.  Within  the  framework  of  the  first  criterion  of  the  economic 
life,  the  percentage  (or  numbers)  of  cracks  exceeding  ag  is  obtained  as  a 
function  of  the  service  time.  The  percentage  of  cracks  exceeding  ae  for 
both  cases,  with  or  without  a  scheduled  maintenance,  is  schematically 
shown  in  Fig.  (2).  For  the  second  criterion  of  the  economic  life,  the 
average  cost  of  maintenance,  including  the  cost  of  inspection  and  repair, 
has  been  formulated  as  a  function  of  service  time,  thus  permitting  a 
determination  of  the  economic  life. 

The  analytical  methodology  is  demonstrated  herein  by  a  numerical 
example.  A  durability  critical  component  of  a  fighter  aircraft  under 
scheduled  inspection  and  repair  maintenance  procedures  has  been  considered 
to  demonstrate  the  methodology  developed  herein.  The  number  of  cracks  in 
the  component  exceeding  the  economic  repair  crack  size  with  any  probabil¬ 
ity  and  confidence  levels  is  obtained  as  a  function  of  service  time.  It 
is  shown  numerically  that  the  number  of  cracks  exceeding  the  economic 
repair  crack  size  (with  any  probability  and  confidence  levels)  increases 
rapidly  after  a  certain  service  time,  thus  determining  the  possible  eco¬ 
nomic  life  of  the  durability  critical  component.  While  the  inspection 
and  repair  maintenance  procedures  have  a  significant  impact  on  the  safety 
and  reliability  of  aircraft  structures  [e.g.,  11-17],  its  effect  on  the 
economic  life  of  a  component  is  shown  to  be  limited. 


1ST  SERVICE  2ND  SERVICE  3RD  SERVICE  4TH  SERVICE 

INTERVAL  ^  INTERVAL  ^  INTERVAL  ^  INTERVAL  ^ 


Figure  1  Inspection  and  Repair  Maintenance  Schedule 
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CRACK  PERFORMANCE 
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Figure  2  Durability  Analysis  Elements 


SECTION  II 


TECHNICAL  OVERVIEW 


An  important  input  parameter  to  the  durability  analysis  is  the 
initial  fatigue  quality  of  the  critical  parts  of  the  aircraft  struc¬ 
tures.  For  engineering  analysis  and  design  purposes,  the  initial 
fatigue  quality  has  been  characterized  by  either  the  time  to  crack  ini¬ 
tiation  (TTCI)  or  the  equivalent  initial  flaw  size  (EIFS) ,  and  attempts 
have  been  made  in  the  literature  [6,  18-35]  to  characterize  the  statis¬ 
tical  distributions  of  both  using  available  laboratory  or  field  data. 

The  time  to  crack  initiation  is  defined  as  the  time  at  which  a  vis¬ 
ible  crack  size  ag,  say  0.03",  occurs.  In  the  present  development  of 
the  durability  analysis  methodology,  the  input  initial  quality  is  chosen 
to  be  the  time  to  crack  initiation  (TTCI) .  The  reasons  are: (i)  the 
crack  size  ag  at  crack  initiation  is  physically  observable,  and  (ii) 
test  data  can  be  obtained  from  the  coupon  specimens  and  the  full-scale 
components,  including  the  test  data  of  fractography . 

The  time  to  crack  initiation,  however,  is  a  function  of  so  many 
variables  that  if  any  one  variable  changes,  such  as  loading  spectra  or 
stress  level,  test  data  may  have  to  be  generated  again,  since  the  func¬ 
tional  relation  between  TTCI  and  other  influencing  variables  is  not 
available  to  date.  This  situation  becomes  particularly  critical  in  the 
present  durability  analysis  where  the  crack  growth  damage  at  each  loca¬ 
tion  of  the  entire  durability  component  has  to  be  estimated.  Since  the 
maximum  stress  level  of  the  loading  spectra  may  vary  from  one  location 
to  another,  it  may  not  be  economically  practical  to  conduct  laboratory 
tests  for  TTCI  for  all  maximum  stress  levels  which  may  occur  at  all  the 
fastener  holes. 

As  a  result,  the  input  data  of  TTCI,  which  is  statistically  charac¬ 
terized  by  the  three-parameter  Weibull  distribution  [20,  21], is  converted 
(transformed)  through  backward  extrapolation  herein  by  use  of  a  crack 
propagation  law  to  determine  the  statistical  distribution  of  the  "equi¬ 
valent"  initial  flaw  size  [20,21,36]  expediently.  Such  an  approach  is 
motivated  by  the  expedient  necessity  stated  above.  Various  questions 
associated  with  such  an  approach  will  be  justified  later,  except  to 
emphasize  that  the  equivalent  initial  flaw  sizes  thus  obtained  for  one 
type  of  particular  aircraft  cannot  necessarily  be  used  for  other  types 
of  aircraft  with  completely  different  mission  characteristics.  Note 
that  the  statistical  distribution  of  the  equivalent  initial  flaw  size 
derived  from  the  three-parameter  Weibull  distribution  of  TTCI  is  consis¬ 
tent  with  the  physical  fatigue  wear-out  process  of  metals  [20,21,36]  . 
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It  is,  therefore,  very  important  to  emphasize  that  the  durability  analysis 
methodology  developed  herein  does  not  distinguish  the  time  to  crack  ini¬ 
tiation  approach  from  that  of  the  equivalent  initial  flaw  size  approach. 
This  is  because  both  approaches  are  interrelated  .  Again,  it  should  be 
noted  that  one  set  of  the  equivalent  initial  flaw  size  data  obtained 
from  one  type  of  aircraft  cannot  necessarily  be  used  for  other  types  of 
aircraft  with  different  mission  characteristics. 

The  entire  population  of  the  equivalent  initial  flaw  size,  as  speci¬ 
fied  by  its  statistical  distribution  derived  from  TTCI  data, is  then  sub¬ 
jected  to  service  loading  spectra  consisting  of  possibly  various  different 
missions.  The  crack  growth  damage  accumulation  under  spectrum  loading  is 
computed  using  a  general  crack  growth  master  curve  which  can  be  obtained 
either  by  a  general  computer  code  or  by  any  crack  growth  model  or  experi¬ 
mental  results,  including  those  of  f ractography .  It  should  be  men¬ 
tioned  that  the  crack  growth  damage  calculation  approach  recently  proposed 
by  Gallagher  [37-39]  can  also  be  used  conveniently. 

Thus,  the  statistical  distribution  of  crack  sizes  at  any  service 
time  T  can  be  obtained  from  that  of  the  equivalent  initial  flaw  size 
distribution  through  the  transformation  of  the  crack  growth  damage 
accumulation  master  curve.  Then,  the  percentage  of  cracks  exceeding 
the  economic  repair  crack  size  ae  with  certain  probability  and  con¬ 
fidence  levels  can  be  computed  as  shown  in  Fig.  (2). 

A  maintenance  program  can  be  implemented  in  service  where  inspections 
and  repairs  are  performed  at  scheduled  intervals.  The  statistical  uncer¬ 
tainties  of  the  NDI  procedure  in  crack  detection  has  been  accounted  for. 
When  a  cracked  detail  is  detected  and  repaired,  its  fatigue  strength  is 
assumed  to  be  renewed  in  the  sense  that  its  crack  size  distribution  is 
identical  to  that  of  the  equivalent  initial  flaw  size  prior  to  service, 
and  its  crack  propagation  characteristic  remains  unchanged,  i.e.,  neglec¬ 
ting  the  possible  effects  of  cold-work  and  others.  The  crack  growth 
damage  accumulation  of  the  repaired  details,  referred  to  as  the  renewal 
details,  has  also  been  taken  into  account. 

During  inspection,  however,  some  cracked  details  may  not  be  detected 
depending  on  the  resolution  capability  of  a  particular  NDI  procedure 
employed.  When  a  crack  is  missed  during  inspection,  it  continues  to 
grow  until  the  next  inspection,  at  which  time  it  will  have  a  higher  prob¬ 
ability  of  being  detected.  However,  the  probability  of  failure  of  the 
component  increases.  The  reliability  of  the  entire  durability  critical 
component  in  any  service  interval  has  been  formulated  analytically. 

As  a  result,  the  statistical  distribution  of  crack  sizes  is  shifted 
and  changed  continuously  in  service  due  to  crack  propagation.  It  is  also 
subjected  to  truncation  and  modification  after  each  inspection  and  repair. 
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The  statistical  distribution  of  crack  sizes  at  any  service  time  is 
derived,  and  the  number  of  cracked  details  exceeding  any  crack  size 
with  certain  probability  and  confidence  levels  are  obtained  as  a  func¬ 
tion  of  service  time. 

Finally,  the  average  number  of  cracks  detected  and  repaired  during 
each  maintenance  is  obtained  and  the  average  cost  of  maintenance  is 
formulated  as  a  function  of  the  service  life. 
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SECTION  III 

STATISTICAL  CHARACTERIZATION  OF  INITIAL  FATIGUE  QUALITY 


One  of  the  most  important  factors  affecting  the  durability  and 
fatigue  life  of  a  structural  component  is  the  initial  fatigue  quality. 
The  initial  fatigue  quality  of  a  component  depends  on  such  variables  as 
material  properties,  drilling  techniques,  welding  qualities,  manufac¬ 
turing  processes,  assembling  of  structures,  etc.  To  characterize  the 
initial  fatigue  quality  as  a  function  of  various  influencing  variables 
is  a  major  undertaking.  Nevertheless,  for  analysis  and  design  purposes, 
the  initial  fatigue  quality  can  be  characterized  by  either  the  time  to 
crack  initiation  (TTCI)  or  the  equivalent  initial  flaw  size  (EIFS) . 


The  time  to  fatigue  crack  initiation  can  be  defined  as  the  time  at 
which  a  visible  crack  size  a q,  say  0.03",  occurs.  For  practical  engi¬ 
neering  purposes,  attempts  have  been  made  to  characterize  the  statis¬ 
tical  distribution  of  time  to  fatigue  crack  initiation  for  applications 
to  fatigue  analysis  and  design  of  aircraft  structures  [22-34] . 

The  advantages  of  the  time  to  crack  initiation  (TTCI)  approach  are 
(i)  the  crack  size  at  crack  initiation  is  physically  observable,  (ii) 
test  data  can  be  obtained  from  coupon  specimens  and  full  scale  com¬ 
ponents  directly,  and(iii)  linear  fracture  mechanics  for  fatigue  crack 
propagation  can  be  applied  after  crack  initiation.  Hence,  the  fatigue 
process  of  metals  has  been  traditionally  divided  into  three  stages: 

(a)  fatigue  crack  initiation,  (b)  fatigue  crack  propagation  (stable 
crack  propagation),  and  (c)  fatigue  fracture  (unstable  crack  propagation) 

The  disadvantage,  however,  is  that  it  is  a  function  of  several 
variables,  such  as  stress  level,  loading  spectra,  manufacturing  processes 
etc.  Thus,  if  any  one  of  the  influencing  variables  changes,  e.g.,  stress 
level  or  loading  spectra,  test  data  may  have  to  be  generated  again,  since 
the  functional  relation  between  TTCI  and  these  influencing  variables  has 
not  been  established  to  date.  This  situation  becomes  especially  criti¬ 
cal  in  the  durability  analysis  of  aircraft  structures  where  the  crack 
growth  damage  at  each  location  of  the  critical  component  has  to  be  esti¬ 
mated.  Since  the  maximum  stress  level  of  the  loading  spectrum  may  vary 
from  one  location  to  another,  it  is  economically  impractical  to  conduct 
laboratory  tests  for  TTCI  for  each  maximum  stress  level. 
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Another  characterization  of  the  initial  fatigue  quality  is  the 
initial  flaw  size  existing  in  the  structure  prior  to  service  [6,18,19,35]. 
Detectable  flaws  or  defects  do  exist  in  the  components  prior  to  service 
[16,10] .  However,  most  of  the  "initial  flaws"  existing  in  the  structure 
are  not  detectable  by  any  NDI  technique.  Thus,  the  initial  flaw  size  is 
usually  obtained  from  the  results  of  laboratory  fatigue  tests.  When 
detectable  flaws  are  observed  during  fatigue  tests,  they  are  extrapolated 
backwards  using  a  suitable  crack  propagation  law  to  estimate  their  "equi¬ 
valent"  initial  flaw  size  (EIFS)  [6,18,19,35].  The  backward  extrapola¬ 
tion  has  to  be  validated  using  f ractographic  results.  This  approach  implies 
that  the  entire  fatigue  process  is  essentially  subcritical  flaw  growth. 
Test  data  of  the  equivalent  initial  flaw  size  thus  obtained  have  been 
available,  and  their  distributions  have  been  characterized  by  the  four- 
parameter  Johnson  distribution  [6,10,18,19,35]. 

The  advantage  of  this  approach  is  that,  once  the  initial  flaw  size 
distribution  is  established,  the  fatigue  crack  propagation  under  various 
service  loading  spectra  can  be  predicted  analytically  without  further 
experimental  tests.  Furthermore,  the  equivalent  initial  flaw  size 
approach  can  be  used  to  study  the  effect  of  various  factors,  such  as 
drilling  techniques,  manufacturing  processes,  etc.,  on  the  fastener  hole 
quality  [35] . 

The  input  initial  fatigue  quality  to  the  durability  analysis  method¬ 
ology  developed  herein  is  the  time  to  crack  initiation,  because  of  the 
advantages  stated  above  and  because  of  the  familiarity  of  practicing 
engineers.  TTCI  data  should  be  obtained  from  the  test  results  under 
design  loading  spectra.  Since  TTCI  is  a  function  of  several  variables, 
in  particular  the  stress  level  in  the  loading  spectra,  and  since  the 
stress  level  varies  from  one  region  of  the  durability  critical  component 
to  another,  it  may  be  economically  impractical  to  generate  test  results 
of  TTCI  for  each  stress  level.  This  is  critical  to  the  durability  anal¬ 
ysis  in  which  the  entire  crack  population  in  the  whole  component  should 
be  taken  into  account. 

As  a  result,  the  input  TTCI  data  under  the  design  loading  spectra 
with  one  maximum  stress  level  is  used  in  the  analysis.  Then,  the  TTCI 
data  is  proposed  to  be  fitted  by  the  three-parameter  Weibull  distribution 
[20,21],  The  statistical  distribution  of  the  equivalent  initial  flaw 
size  is  derived  from  the  three-parameter  Weibull  distribution  of  TTCI 
and  the  durability  analysis  is  started  from  the  statistical  distribu¬ 
tion  of  the  EIFS.  This  expedient  approach  simplifies  significantly  the 
durability  analysis  procedures  [20]  as  will  be  discussed  later  along 
with  the  justification  of  such  an  approach. 
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Since  the  Weibull  distribution  for  TTCI  represents  the  wear-out 
fatigue  process,  the  derived  new  distribution  for  the  equivalent  initial 
flaw  size  bears  the  physical  notion  of  the  wear-out  fatigue  process. 
Methods  of  estimating  the  distribution  parameters  and  confidence  levels 
will  be  described.  Finally,  test  results  of  TTCI  and  EIFS  are  presented 
to  demonstrate  the  validity  of  the  proposed  distribution.  It  will  be 
shown  that  the  derived  distribution  fits  the  EIFS  data  better  than  any 
other  distribution  function. 


3.1  Time  to  Crack  Initiation  (TTCI) 

Traditionally,  the  lognormal  distribution  has  been  used  to  describe 
the  statistical  distribution  of  the  time  to  fatigue  crack  initiation 
(TTCI)  [22-23] .  However,  the  failure  rate  (or  risk  function)  associated 
with  the  lognormal  distribution  initially  increases  and  then  eventually 
decreases.  This  is  inconsistent  with  the  fatigue  wear-out  model.  From 
a  physical  standpoint,  it  is  generally  agreed  that  metal  fatigue  is  a 
wear-out  process.  For  example,  the  longer  a  specimen  has  survived 
fatigue  testing,  the  probability  of  specimen  failure  increases.  Thus, 
the  fatigue  life  failure  rate  should  increase  monotonically .  If  the 
failure  rate  of  the  fatigue  life  is  a  positive  power  function  of  the  ser¬ 
vice  time,  then  it  can  be  derived  mathematically  that  the  distribution 
function  of  the  time  to  crack  initiation  is  Weibull,  i.e., 

a 

FT(t)  =  P[T<t]  =  1  -  e"(t/6)  ,  t  >  0  (1) 

in  which  T  is  a  random  variable  indicating  the  time  to  crack  initiation, 
a  is  the  shape  parameter  and  3  is  the  scale  parameter.  In  Eq.  (1) , 

F^(t)  =  P[T<t]  is  the  probability  that  the  time  to  crack  initiation  T  is 
smaller  than  or  equal  to  a  value  t. 

The  Weibull  distribution  given  by  Eq.  (1)  has  been  used  to  describe 
the  time  to  crack  initiation  [24-34] ,  and  methods  of  estimating  the  param¬ 
eters  a  and  3  from  test  data  are  available  [26-27]  . 

Observations  of  extensive  specimen  test  data  indicate  that  a  is 
fairly  constant  for  a  particular  material  and  it  is  not  sensitive  to 
specimen  geometry,  testing  method,  specimen  size,  etc.  Compilation  of 
coupon  test  data  indicates  that  a  =  4.0  may  be  appropriate  for  aluminum 
[27] .  Moreover,  attempts  have  been  made  to  apply  the  Weibull  distribu¬ 
tion  for  the  time  to  crack  initiation  to  available  service  data  for  var¬ 
ious  types  of  aircraft,  for  instance  C-130,  C-141,  F4,  etc.  [30-34]. 

An  improvement  in  predicting  the  time  to  service  crack  initiation  has 
also  been  made  by  taking  into  account  the  statistical  variability  of 
service  loads  [34] . 
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The  lower  bound  of  the  two-parameter  Weibull  distribution  given 
by  Eq.  (1)  is  zero.  The  lower  bound,  in  effect,  depends  on  the  defi¬ 
nition  of  the  crack  size  ag  at  crack  initiation.  As  a  result,  the  fol¬ 
lowing  three-parameter  Weibull  distribution  should  be  used  for  gener¬ 
ality  [20]  , 


a 

Fx(t)  =  P[T<t]  =  1  -  exp  {-  (~£)  }  ;  t>£  (2) 

in  which  £  is  the  lower  bound  of  TTCI .  The  lower  bound  £  increases  as 
the  defined  crack  size  ag  at  crack  initiation  increases. 

3.2  Fatigue  Crack  Propagation  in  Small  Crack  Size  Region 

As  mentioned  previously,  the  equivalent  initial  flaw  sizes (EIFS) 
are not  observable  by  any  NDI  technique.  They  are  obtained  from  detec¬ 
table  crack  sizes  observed  during  fatigue  tests  by  backward  extrapola¬ 
tion.  Such  an  extrapolation  must  be  verified  using  fractographic 
results.  A  significant  contribution  to  the  prediction  of  crack  growth 
damage  under  complex  flight  loading  has  been  made  recently  by  Gallagher 
who  indicates  that  the  crack  growth  rate  under  flight-by-flight  loading 
spectra  can  be  written  as  follows  [37-39]  . 

ft  -G(5WA  (3) 

in  which  t  represents  the  flight  (or  flight  hour  or  miniblock  of  spectra) , 
G  and  X  are  parameters  depending  on  such  factors  as  loading  spectra, 
material  and  structural  properties,  etc. 

In  Eq.  (3),  K  is  the  root  mean  square  of  the  maximum  stress 
.  .  .  ^  ^  max 

intensity  factor. 


where  q  is  the  mean  square  of  the  maximum  stress  in  the  flight  spectra 

and  ^(a)  can  be  expressed  in  terms  of  the  stress  intensity  factor  coeffi¬ 
cient  3(a) ,  as 
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4*  (a)  =  3(a)v/TTa 


(5) 


The  stress  intensity  factor  coefficient  3(a)  depends  on  the  crack 
size,  the  specimen  width,  the  crack  geometry  (e.g.,  through -the -thick¬ 
ness  crack,  part-through  crack,  corner  crack),  etc.  In  general, 

ip(a)  can  be  approximated  by  a  power  series  of  the  crack  size  a  as 

oo  m 

ifj(a)  =  £  n.  a  (6) 

i=l  1 

In  the  region  of  small  crack  size,  i.e.,  smaller  than  the  crack 
size  aQ  at  crack  initiation,  a  good  approximation  can  be  obtained  by 
taking  the  first  term  only,  i.e.. 


<Ka) 


(7) 


Substituting  Eqs.  (4)  and  (7)  into  Eq.  (3)  ,  one  obtains 


=  Qab(t)  (8) 

The  advantage  of  the  miniblock  approach  proposed  by  Gallagher 
[37-39]  is  that  the  crack  growth  damage  can  be  computed  by  performing 
the  flight-by-flight  integration  rather  than  the  cycle-by-cycle  inte¬ 
gration,  thus  reducing  a  significant  amount  of  computer  time. 

While  the  validity  of  the  miniblock  approach  has  been  verified 
[37-39],  Eq.  (8)  needs  to  be  verified  using  physical  observations. 

If  plots  are  obtained  of  log  da/dt  versus  log  a  for  all  specimens 
tested  in  a  data  set,  then  a  curve  fitting  technique  can  be  used  to  find 
the  best-fitting  straight  line  through  the  data  points.  A  least  squares 
technique  is  sufficient  and  by  minimizing  the  least  squares  error, 
"composite"  values  of  b  and  Q,  the  slope  and  intercept  of  the  line  ,  respec¬ 
tively,  can  be  determined  for  the  data  set  as  a  whole.  Typical  plots 
of  test  data  and  Eq.  (.8)  demonstrate  good  agreement  as  shown  in  Figs. 

(3)  and  (4).  The  relation  given  by  Eq.  (8)  eventually  becomes  invalid 
when  the  crack  size  becomes  large  as  will  be  discussed  later. 


13 


•  •  • 
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Figure  4  Log  Crack  Length  Versus  Log  da/dt  for  XWPF  (37)  Data  Set 


By  back  extrapolation  from  the  observed  TTCI  values  using  Eq.  (8) 
and  these  composite  Q  &  b  values,  the  EIFS  distribution  can  be  found. 
Although  this  distribution  cannot  be  verified  by  direct  physical  observ¬ 
ation  ,  when  grown  out  to  some  time  t  again  using  Eq.  (8)  ,  the  predicted 

distribution  of  crack  sizes  at  this  time  t  compares  very  closely  with  the 
observed  distribution.  This  suggests  that  the  crack  growth  model  is 
valid  and  that  the  composite  Q  and  b  values  obtained  represent  an  aver¬ 
age  growing  crack  sufficiently  well.  These  findings  are  detailed  in 
Section  3.5. 


3.3  Statistical  Distribution  of  EIFS 


Let  a(T)  be  the  crack  size  at  T  flight  hours,  where  T  is  the  time 
to  crack  initiation,  and  a(0)  be  the  initial  crack  size  (at  t=0)  prior 
to  service.  It  is  obvious  that  a(T)=aQ  is  the  crack  size  at  crack  initia¬ 
tion.  Integration  of  Eq.  (8)  from  t=0  to  t=T  yields  the  relation 
between  the  initial  flaw  size  a(0)  and  Sq, 


EIFS  = 


a(0) 


(1+Sq  CQT)1/C 


(9) 


in  which 

c  =  b  -  1 


(10) 


The  statistical  distribution  of  time  to  crack  initiation  T,  is 
given  by  Eq.  (2)  .  Thus,  the  statistical  distribution  of  the  initial 

flaw  size,  a(0),  can  be  derived  from  Eq.  (2)  through  the  transforma¬ 
tion  of  Eq.  (9)  as  follows: 


(ID 


Substituting  Eq.  (2)  into  Eq.  (11)  ,  one 

function  F  (x)  of  the  initial  crack  size  as 
a(0) 


obtains 

follows 


the  distribution 
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a(0) 


— c  —  c 

a) 

exp  j- 

X  -a  -cQe 

cQ8 

J 

;  ) 


=  1; 


x>x 
—  u 


(12) 


where  x  is  the  upper  bound  of  the  initial  crack  size. 

/  -c  ^  \ -1/c 
*u  ’  (a0  +cQe  ) 

Such  a  transformation  is  schematically  shown  in  Fig.  (5). 


(13) 


The  EIFS  cumulative  distribution,  F  (x)}Eq.  (12),  cannot  be 

a(.u) 

verified  directly  from  f ractographic  results.  However,  the  cumulative 
distribution,  F  (x) ,  of  the  crack  size  a(t)  at  any  service  time  t 

where  crack  sizes  are  observable  from  fractographic  results  can  be 

derived  from  F  ...  (x).  A  verification  of  F  .  (x)  then  gives  an  indirect 

a  (.  t) 

verification  of  F  ...  (x) .  F  .  .  (x)  is  verified  using  actual  fractog- 
a(0)  a(t) 

raphy  in  Section  3.5. 


The  cumulative  distribution  of  crack  sizes  at  any  service  time  t, 

F  .  . (x) ,  is  described  by  Eqs.  (14)  and  (15)  (See  Appendix  for 
a(t) 

derivation) , 


(14) 


where, 


yx(t) 


_ X _ 

[l+xCcQt] 


(15) 


3.4  Estimation  of  Parameters  and  Confidence  Level 

The  three  parameters,  a,  8>  and  e,  appearing  in  Eq.  (2)  should  be 
determined  from  test  results  of  time  to  crack  initiation.  The  curve 
fitting  procedure  along  with  the  least  square  method  can  be  employed  to 
estimate  a,  8,  and  e  as  follows. 
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CRACK  SIZE 


Figure  5  Elements  of  Initial  Fatigue  Quality  Model 
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First,  a  value  is  assumed  for  the  location  parameter  £.  Then, 
Eq.  (2)  is  transformed  into  the  following 


£n[-£n(l-F  (t))]  =  a£n(t-e)-a£n£ 
T 

Let 

Z  =  £n[-£n(l-Fx(t))] 

Y  =  £n(t-e) 

U  =  -cdn3 


(16) 


(17) 


Then,  Eq.  (16)  reduces  to  a  straight  line  in  the  (Y,Z)  plane, 
referred  to  as  the  probability  paper, 

Z  =  aY  +  U  (18) 

Test  data  of  time  to  crack  initiation  can  be  plotted  on  the  proba¬ 
bility  paper,  and  a  and  3  are  estimated  using  the  least  square  fit. 

Another  value  of  £  is  then  assumed  and  the  same  procedure  is  repeated 
to  estimate  another  set  of  a  and  3-  The  final  values  of  a,  3,  and  £  are 
chosen  to  minimize  the  corresponding  mean  square  error  between  the  test 
data  and  the  fitted  distribution. 

For  given  a  and  e,  the  y  confidence  level  for  3,  denoted  by  3y>  is 
given  by  [27] 


in  which  n  is  the  total  number  of  test  data  points  and  K^,(2n)  is  the  y- 
fractile  of  the  chi-square  variate  with  2n  degrees  of  freedom.  Eq.  (19) 
would  be  exact  If  the  estimate  of  3  were  obtained  using  the  maximum- 
likelihood  method. 


3.5  Test  Results  And  Correlation 

Test  results  from  the  "Fastener  Hole  Quality"  program  [35]  will  be 
used  to  verify  the  EIFS  distribution,  Eq.  (12)  •  Test  results  for  two 

types  of  Aluminum  (7475-T7351)  specimens  subjected  to  a  F-16  fighter  400-hour 
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block  spectrum  with  a  maximum  gross  stress  of  34  ksi  are  used:  (1)  WPF 
(no  load  transfer)  and  (2)  XWPF  (15%  load  transfer) .  Specimen  details 
are  shown  in  Fig.  (6)  .  Observed  TTCI's  at  a  crack  size  of  0.03"  along 
with  crack  sizes  at  two  different  times  for  WPF  and  XWPF  data  sets 
are  presented  in  Tables  1  and  2,  respectively. 

The  WPF  (38)  data  set  (Table  1)  includes  five  specimens  with  only 
one  fractographic  observation/specimen  (i.e.,  one  crack  size  and  one 
TTCI) .  Crack  sizes  and  TTCI’s  were  computed  for  these  five  WPF  specimens 
(Table  1)  using  Eq.  (8)  and  the  crack  growth  model  constants,  Q  and 
b  (Table  3),  based  on  the  WPF  (33)  data  set.  The  WPF  (33)  data  set 
includes  all  specimens  in  Table  1  with  observed  TTCI's  and  applicable 
crack  sizes.  Both  data  sets,  WPF  (38)  and  WPF  (33),  are  used  later. 

The  EIFS  cumulative  distribution  will  be  validated  as  follows: 

First,  it  will  be  shown  that  the  three-parameter  Weibull  distribution 
fits  the  observed  TTCI  results  for  a^=0.03".  Second,  the  theoretical 
cumulative  EIFS  distribution  will  be  compared  against  ranked  EIFS  pre¬ 
dictions  (t=0).  Third,  the  predicted  cumulative  crack  size  distributions 
for  two  different  times  and  two  different  specimens  will  be  compared 
against  the  corresponding  ranked  crack  sizes  observed. 

Weibull  best  fit  plots  for  WPF  (38),  WPF(33),  and  XWPF  data  sets 
are  shown  in  Figs.  (7)  through  (9)  ,  respectively.  The  three-parameter 
Weibull  distribution  fits  the  observed  TTCI  results  well.  Parameter 
values  used  are  summarized  in  Table  3. 


The  values  of  a,  8  and  e  from  the  test  results  of  TTCI  (Tables  1 
and  2)  and  the  values  of  b  and  Q  based  on  crack  growth  fractography  data 
(Figs.  (3)  and  (4).),  can  now  be  used  to  examine  the  correlation  between 
the  observed  and  the  theoretically  predicted  results. 

EIFS  values  are  computed  for  three  data  sets,  WPF  (38),  WPF  (33) 
and  XWPF,  using  the  TTCI  results  of  Tables  1  and  2  and  Eq.  (9)  .  Results 

are  ranked  and  plotted  in  Figs.  (lO)through  (12)  as  circles.  Predicted 
cumulative  distribution  curves,  represented  by  Eq.  (12)  ,  are  plotted 

in  Figs. (10)  through  (12)  .  Furthermore,  observed  crack  sizes  at  two 
times,  presented  in  Tables  1  and  2,  are  also  ranked  and  plotted  in  Figs. 
(10)  through  (12).  and  compared  with  the  predicted  cumulative  distribution 
represented  by  Eq.  (14)  for  WPF  (38),  WPF  (33),  and  XWPF  data  sets.  The 
correlation  between  the  observed  and  the  predicted  results  is  excellent. 
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NO  LOAD  TRANSFER  SPECIMEN 


d  -  0.250" 


15%  LOAD  TRANSFER  SPECIMEN 


d  »  0.250" 


Figure  6  Test  Specimen  Details 
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Table  1  Observed  Test  Results  For  WPF  Data  Set 


SPECIMEN 

NO. 

0 

WPF  (3 eP 

Hi 

a(t  =  9200  Hrs) 

(Inch) 

IB1 

WPF  -11 

6563 

0.0617 

0.2570 

-  7 

6312 

0.0287 

0.3557 

-16 

10629 

0.0229 

0.0582 

-21 

10892 

0.0205 

0,0654 

-  9 

11357 

0.0185 

0.0570 

-  6 

11637 

0.0142 

0.0687 

-  8 

12053 

0.0143 

0.0548 

-18 

13379 

0.0118 

0.0387 

-23 

13544 

0.01  IS 

0.0378 

-32 

13762 

0.0099 

0.0372 

-14 

13773 

0.0112 

0.0357 

-17* 

(13783) 

(0.0067) 

(0.0392) 

-25 

13339 

0.00B2 

0.0360 

-35 

14098 

0.0072 

0.0354 

-13 

14123 

0.0082 

0.0348 

-12 

14149 

0.0064 

0.0361 

-19 

14252 

0.0073 

0.0332 

-22 

14350 

0.0082 

0.0329 

-10 

0.0082 

0.0322 

-34 

0.0080 

0.0322 

-29 

iBlIi 

0.0056 

0.0259 

-24 

15600 

0.0049 

0.0252 

-33 

15798 

0.0038 

0.0225 

-20 

16128 

0.0034 

0.0217 

16141 

0.0047 

0.0215 

-31 

17109 

0.005S 

0.0185 

-37 

17134 

0.0033 

0.0172 

-38 

17185 

0.0048 

0.0170 

-38 

17507 

0.0052 

0.0162 

-36 

17620 

0.0032 

0.0151 

-26 

17639 

0.0037 

0.0148 

-30 

18066 

0.0037 

0.0128 

-27 

0.0028 

0.0121 

-41* 

(18696) 

(0.0022) 

(0.0105) 

-28 

0.0031 

0.0099 

-43* 

(20140) 

0.0015 

0.0071 

-15* 

(20449) 

(0.0013) 

(0.0065) 

-42* 

(0.D005) 

(0.0B26) 

NOTES 

(D  Ref.  Fig.  6 
©For  aQ  =  0.03" 

©WPF(38)  =  total  WPF  data  set 

*  Single  observation/specimen  {i.e.,  single  data  point:  a, ,  t-) 
aj  =  crack  size 

tj  =  time 

(rax)  =  Value  derived  from  single  observation  using  Eq.  8  and  results  in 
Table  3  for  WPF(33) 
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Table  2  Observed  Test  Results  For  XWPF  Data  Set 


XWPF  (Ref.  35) 

a(t  -  7200  Hrs) 
(Inch) 

a(t  -  10800  Its) 
(Inch) 

XWPF- 36 

6106 

0.0392 

0.0814 

-38 

6147 

0.0418 

0.0922 

-18 

7457 

0.0264 

0.0943 

-  9 

7545 

0.0248 

0.0940 

-  7 

7721 

0.0211 

0.1098 

-  3 

8425 

0.0186 

0.0910 

-34 

8528 

0.0153 

0.0753 

-31 

8636 

0.0195 

0.0525 

-29 

8908 

0.0206 

0.0429 

-15 

8968 

0.0140 

0.0582 

-  8 

9078 

0.0168 

0.0487 

-22 

9085 

0.0131 

0.0575 

-  1 

9316 

0  0195 

0.0386 

-  2 

9973 

0.0111 

0.0402 

-11 

10253 

0.0067 

0.0363 

-12 

10457 

0.0081 

0.0333 

-16 

10908 

0.0088 

0.0288 

-24 

11045 

0.0120 

0.0281 

-19 

11051 

0.0049 

0.0271 

-30 

11071 

0.0052 

0.0257 

-  5 

11370 

0.0087 

0.0264 

-23 

11493 

0.0042 

0.0201 

-14 

11564 

0.0022 

0.0110 

-28 

11571 

0.0047 

0.0213 

-17 

11698 

0.0063 

0.0223 

-27 

11708 

0.0058 

0.0210 

-37 

11920 

0.0048 

0.0215 

-33 

12051 

0.0071 

0.0203 

-13 

12062 

0.0045 

0.0174 

-  6 

12118 

0.0061 

0.0208 

-21 

12302 

0.0060 

0.0191 

t25 

12389 

0.0070 

0.0193 

-26 

12505 

0.0040 

0.0140 

-35 

12629 

0.0045 

0.0169 

-  4 

12655 

0.0041 

0.0136 

-10 

13339 

0.0030 

0.0100 

-20 

16008 

0.0007 

0.0034 

NOTES:  ®  Ref.  Fig.  6 
©  aQ  *  0.03” 
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Table  3  Model  Parameters  Based  On  Observed  Test  Results 


Parameter 

Specimen 

WPF(38) 

WPF(33) 

XWPF(37) 

b 

0.9703 

0.9703 

0.9620 

c 

-0.0297 

-0.0297 

-0.0380 

Q 

0.0002381 

0.0002381 

0.000309 

Of 

4.863 

4.9174 

5.499 

0 

14,957 

15,936 

11,193 

Py 

14,240 

15,102 

10,694 

Y 

€ 

1,312 

0 

0 
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Figure  7  Weibull  Best  Fit  Plot  for  WPF  (38)  Data  Set 
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Ref  35 
(33  Specimen) 
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Figure  8  Weibull  Best  Fit  Plot  for  WPF  (33)  Data  Set 


In  -In  1-FlW 


2. 


Y=ln  (t— 0 


Figure  9  Weibull  Best  Fit  Plot  for  XWPF  Data  Set 
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50%  CONFIDENCE 
95  %  CONFIDENCE 


Figure  10  Crack  Size  Versus  Cumulative  Distribution 
F  (x)  for  WPF  (38)  Data  Set 


5 


(x)a:°j  Nounamisia  BAiiviniAino 
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Figure  11  Crack  Size  Versus  Cumulative  Distribution 
F  ,^00  for  (33>  Data  Set 


CUMULATIVE  DISTRIBUTION 


Figure  12 


Crack  Size  Versus  Cumulative  Distribution  F  ,  .  (x) 
for  XWPF  Data  Set 


Other  statistical  distributions,  such  as  the  normal,  lognormal,  Wei- 
bull,  Pearson,  and  Johnson, have  been  used  to  fit  the  test  results  [36].  The 
Kolmogorov- Smirnov  test  for  goodness  of  fit  indicates  the  presented 
theoretical  distribution,  Eq.  (12)  ,  referred  to  as  the  Weibull- 

compatible  distribution,  fits  the  test  results  best. 


The  predicted  cumulative  distribution  of  crack  sizes  F  ,  . 

a(t) 


(x) 


tracks 


observed  crack  sizes  very  well  (Figs.  (10)  through  (12) )  .  The  EIFS  dis¬ 
tribution  at  t=0  is  not  only  compatible  with  the  observed  TTCI  results, 
but  it  is  also  statistically  compatible  with  the  observed  crack  size  dis¬ 
tributions,  Although  EIFS  cracks  are  not  observable,  the  above  approach 
provides  an  indirect  verification  of  the  proposed  EIFS  distribution, 

Eq-  (12). 


3*6  Conclusions 

A  physically  meaningful  EIFS  distribution  has  been  suggested  herein. 
The  distribution  is  derived  from  the  three-parameter  Weibull  distribution 
that  is  used  for  the  distribution  of  time  to  crack  initiation.  Since 
the  Weibull  distribution  satisfies  the  condition  of  the  fatigue  wear-out 
process  [24-27] ,  the  EIFS  distribution  suggested  herein  bears  the  physi¬ 
cal  meaning  of  the  fatigue  wear-out  process.  It  is  shown  that  the  cor¬ 
relation  between  the  predicted  crack  size  distribution  at  time  t,  extrapo¬ 
lated  from  the  EIFS  distribution,  and  the  observed  crack  sizes  at  time  t 
is  excellent. 
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SECTION  IV 


CRACK  GROWTH  DAMAGE  ACCUMULATION 


After  the  statistical  distribution  of  the  equivalent  initial  flaw  size  is 
established  from  the  test  results  of  TTCI  described  previously,  the  entire 
population  of  the  EIFS  is  subjected  to  crack  propagation  in  service.  Because 
of  the  crack  geometry,  the  effect  of  loading  sequence  (e.g.,  retardation  and 
acceleration),  and  many  other  factors,  the  crack  growth  rate  equation  given 
by  Eq.  (8)  does  not  hold  for  the  entire  region  of  the  crack  size.  As  a 
result,  crack  growth  predictions  must  be  carried  out  numerically  by  use  of  a 
cycle-by-cycle  integration  approach  using  a  general  computer  program.  This 
program  must  account  for  all  the  variables  influencing  the  crack  growth  behav¬ 
ior  such  as  crack  geometry,  load  sequence  effect,  and  many  others.  The  crack 
size  a(t^)  at  flight  hours  can  be  expressed  in  terms  of  a(t^),  where  t^-t2 
as  follows,  [e.g..  Ref.  6]  : 


a(t2)  =  a (tx)  +  l  Aa(tj) 


(20) 


in  which  Aa(t.)  is  the  crack  growth  increment  per  flight 

t-  <r  J 

1  j  2' 


hour  at  t .  where 
J 


Thus,  the  crack  growth  damage  a(t)  as  a  function  of  the  service  time  t 
is  determined  from  the  general  computer  program  starting  from  a  crack  size  much 
smaller  than  the  EIFS.  Such  an  analytical  crack  growth  curve  is  referred  to  as 
the  "master  curve". 


Two  master  curves  for  the  fastener  holes  of  747 5-T7 351  aluminum  specimens  with 
no  load  transfer  and  Winslow  drilled  with  proper  drilling  techniques  under  a 
fighter  spectrum  are  presented  in  Fig.  (13).  Curves  1  and  2  are  associated 
with  different  maximum  stress  levels  amax  in  the  fighter  spectrum  as  shown  in 
the  figure.  Furthermore,  Fig.  (14)  indicates  two  master  curves  for  the  same 
fastener  hole  conditions  and  the  same  maximum  stress  levels  except  with  a  15% 
load  transfer.  Both  Figs.  (13)  and  (14)  are  established  using  the  genera] 
crack  growth  damage  accumulation  computer  code  developed  by  General  Dynamics. 

The  crack  growth  rates  corresponding  to  Figs.  (13)  and  (14)  are  presented 
in  Figs.  (15)  and  (16)  respectively.  It  is  observed  from  these  figures  that  in  the  small 
crack  size  region  the  crack  growth  rate  vs.  the  crack  size  in  log  scale  is 
practically  a  straight  line,  indicating  the  validity  of  Eq.(8). 
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No  Load  Transfer 


(HDNI)  3ZIS  HYH 
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Figure  13  Flaw  Size  Versus  Flight  Hours  (No  Load  Transfer) 
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Figure  14  Flaw  Size  Versus  Flight  Hours  (15%  Load  Transfer) 


FLAW  GROWTH  RATE  (INCH  PER  FLIGHT  HOUR) 


Figure  15  Flaw  Growth  Rate  Versus  Flaw  Size  (No  Load  Transfer) 
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FLAW  GROWTH  RATE  .(INCH  PER,  FLIGHT  HOUR) 


Figure  16  Flaw  Growth  Rate  Versus  Flaw  Size  (15%  Load  Transfer) 
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It  follows  from  the  particular  case  of  Eq.  8  and  Eqs . (A-l) and (A-2)  f rom  the 
Appendix  that  the  relationship  between  the  crack  sizes  a(t^)  and  a(t2)  at  and 
t  2  does  not  depend  on  t-^  and  t2  but  only  on  t2~t-^.  This  relationship  is  not 
valid  for  the  general  case  within  one  flight  of  a  particular  mission,  because 
of  the  crack  size  and  the  load  sequence  effect.  Thus,  within  one  flight  of  a 
particular  mission,  the  relationship  between  a(t^)  and  a(t2)  represented  by 
Eq.  (20)  depends  on  both  t^  and  t2*  Since,  however,  the  design  loading 
spectra  in  one  lifetime  consists  of  many  repeated  missions  and  flights,  it  is 
reasonable  to  assume  that  the  relationship  between  a(t^)  and  a(t2)  depends  only 
on  the  time  difference  t2~t^  of  the  service  time.  Such  an  expedient  approxima¬ 
tion  appears  to  be  acceptable  and  it  simplifies  the  computational  procedures 
of  the  durability  analysis  tremendously.  As  a  result,  only  a  crack  growth 
master  curve  for  each  maximum  stress  level  in  the  loading  spectra  is  sufficient 
for  the  estimation  of  the  crack  growth  damage  accumulation  in  the  durability 
analysis. 

Thus  for  the  purpose  of  mathematical  derivation,  the  analytical  master  curve 
a(t)  in  the  ith  stress  region  can  be  symbolically  represented  by  the  following 
equation  in  terms  of  t^-t^ 

a(tx)  =W[a(t2),  t2-tJL]  (21) 


in  which  W  is  a  general  function  representing  the  master  curve.  Eq.  (21)  indi¬ 
cates  that  for  a  given  value  of  crack  size  a(t^)  at  t2’  °ne  Can  c*eterm:’-ne  t^ie  crac^ 
size  a(t^)  at  t^  from  the  master  curve  as  depicted  in  Fig.  (17).  Eq.  (21)  will 
be  used  in  the  mathematical  development  of  the  durability  analysis  for  predicting 
the  crack  growth  damage  accumulation  of  the  crack  population.  It  will  serve, 
mathematically,  as  the  transfer  function  for  transferring  the  statistical  dis¬ 
tributions  of  the  crack  population  from  one  service  time  to  another,  as  will  be 
described  later. 

It  should  be  emphasized  that  in  the  present  analysis,  any  general  crack  growth 
damage  accumulation  package  and  method  [6,7,37-47]  can  be  used  to  obtain  the  master 
curve.  It  appears  that  the  miniblock  approach  proposed  by  Gallagher  [37-39]  using 
the  flight-by-flight  integration  technique  is  most  efficient  for  the  present  pur¬ 
pose.  Furthermore,  the  master  curve  varies  from  one  stress  region  to  another  in 
the  entire  durability  critical  component.  Hence,  an  appropriate  master  curve 
should  be  used  for  a  particular  stress  region. 
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SECTION  V 


CRACK  EXCEEDANCE  AND  MAINTENANCE  COST 


A  durability  critical  component  is  divided  into  m  stress  regions.  The  maxi¬ 
mum  stress  level  in  each  stress  region  at  every  location  or  detail  (e.g.,  fastener 
hole)  is  approximately  the  same,  whereas  the  maximum  stress  level  varies  from  one 
stress  region  to  another.  Awing  panel  of  aF-16  fighter  aircraft,  which  is  iden¬ 
tified  as  a  durability  critical  component,  is  shown  in  Fig.  (18),  The  lower 
skin  of  the  wing  panel  is  divided  into  three  stress  regions  where  the  maximum 
applied  stresses  (limit)  to  each  region  under  the  fighter  spectrum  are 
24.3  ksi,  27.0  and  28.3  as  shown  in  Fig.  (19), 

Let  N.  be  the  total  number  of  details  in  the  ith  stress  region.  For  instance, 
the  total  number  of  fastener  holes  in  each  region  of  Fig.  (19)  is,  respectively, 
59,  335  and  1220.  If  the  maximum  stress  level  varies  from  one  location  to  another, 

then  N  .=1. 

l 

In  the  ith  stress  region,  let  N(i,T)  represent  the  total  number  of  details 
having  a  crack  size  exceeding  x^  at  any  service  time  x.  It  is  obvious  that  N(i,x) 
is  a  statistical  (random)  variable,  since  the  initial  flaw  size  is  a  statistical 
variable.  It  is  mentioned  that  each  detail  such  as  the  fastener  hole  may  have 
multiple  cracks.  In  this  report,  the  crack  of  a  detail  refers  to  the  largest  crack, 
and  hence,  each  detail  has  one  crack  (the  largest  one).  As  a  result,  the  term 
"crack"  is  interchangeable  with  "detail". 

It  is  reasonable  to  assume  that  the  crack  growth  at  each  detail  is  not  influ¬ 
enced  by  the  cracks  in  its  neighboring  details  and  that  the  crack  damage  accumula¬ 
tion  of  each  detail  is  statistically  independent.  Then,  the  statistical  distri¬ 
bution  of  N(i,x)  can  be  shown  to  follow  the  Binomial  distribution 

r  -1  /{J  ,\  r  -jN.  "'ll 

P  |n  ( i  #  T )  =  ij  =  pn(i,T)  |JL-p(i,T)J  1  (22) 

in  which  P[N(i,x)  =  n]  denotes  the  probability  that  the  total  number  of  details 
having  a  crack  size  exceeding  x  in  the  ith  stress  region,  N(i,x),  is  equal  to  n, 
and  p(i,x)  is  the  probability  that  a  detail  in  the  ith  stress  region  will  have  a 
crack  size  greater  than  x^  at  the  service  time  x.  In  Eq.  (22) 
nation  of  n  out  of  N^,  i.e., 

HA  N±  1  r(N±+l) 

n  /  nTTNpnJ  f  =  r  (n+l)  r  (N^n+1)" 


■(:■) 


is  the  combi- 


(23) 
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LEGEND:  Approximate  Number  of  Fasteners  Per  Stress  Range 


27.0-28.3  ksi  Limit  (59  Holes) 


24.3-27.0  ksi  Limit  (335  Holes) 
t —  1  Less  Than  21.6-24.3  k<i  I  imit 


Figure  19  Stress  Zoning  for  Durability  Component  Wing-Lower  Skin 
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in  which  T(.)  is  the  gamma  function. 

2 

The  average  value,  N(i,x),  and  the  variance,  a^Ciji),  of  N(i,x)  (the  total 
number  of  details  in  ith  stress  region  having  a  crack  size  greater  than  x^  at 
any  service  time  x)  can  be  obtained  from  Eq.  (22)  as  follows 


E(  i,x) 


=  Z  nP{N(i,x)  =  n] 
n=0 

00  2 

=  Z  n  P[N(i,x>  *=  n] 
n=o 

88  Nip(i,x)  [l-p(i,x)J 


=  N^d^x) 
-  N2(i,x) 


(24) 


(25) 


If  N*  denotes  the  total  number  of  details  in  the  entire  durability  critical 
component  and  L(x)  indicates  the  total  number  of  details  in  the  entire  component 
having  a  crack  size  greater  than  x^  at  the  service  time  x,  then  it  is  obvious  that 


m 


N*  =  Z  N. 

i=l  1 


(26) 


m 


L(x)  =  Z  N(i,x) 
i=l 


(27) 


Thus,  the  average  value,  L  (x) ,  and  the  variance  a  (x) ,  of  L(x)  are  obtained, 
respectively, 


IQ 


L(x)  =  Z  N(i,x) 
i=l 


(28) 


o  m  ^ 

Oj  (x)  =  Z  oN(i,x) 
L  i=l  N 


(29) 


in  which  N(i,x)  and  a^(i,x)  are  given  by  Eqs.  (24)  and  (25)  . 
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The  total  number  of  details  in  the  entire  component,  L(t),  having  a  crack 
size  larger  than  x^  at  the  service  time  T  is  a  statistical  variable.  It  follows 
from  Eq.  (27)  that  L(t)  is  the  sum  of  independent  Binomial  variables,  Eq.  (22), 


When  the  number  of  details  in  each  stress  region,  N.,  is  large,  the  Binomial 
distribution  given  by  Eq.  (22)  can  be  approximated  by  tfte  Normal  distribution, 
in  which  case,  the  statistical  distribution  of  L(t)  is  also  Normal  with  the  mean 
and  the  variance  given  by  Eqs.  (28)  and  (29)  .  The  approximation  by  the  Normal 
distribution  does  simplify  the  computational  effort. 


The  component  may  undergo  a  scheduled  inspection  and  repair  maintenance  pro¬ 
cedure  at  T^,  T^,  T^  ...  in  service  with  service  intervals  T^,  •••  as 

shown  in  Figure  (1).  It  is  assumed  that  the  cracked  details  detected  during 
inspection  will  be  repaired.  Let  g.(x)dx  be  the  probability  of  detecting  (or 
repairing)  a  crack  in  the  size  rang^  (x,x  +  dx)  at  the  jth  inspection.  Then  the 
average  maintenance  cost,  consisting  of  the  cost  of  inspections  and  repairs,  in 
the  service  interval  (0,T^.+^),  denoted  by"C(I),  is 


I 

C(I)  =  C,N*I  +  Z 

1  j=l 


(x)g^(x)dx 


(30) 


in  which  =  cost  of  inspecting  one  detail,  I  =  total  numbers  of  scheduled  inspec¬ 
tions,  C2(x)  =  cost  of  repairing  one  cracked  detail  having  a  crack  size  x,  and 
N*  is  the  total  number  of  details  in  the  entire  component,  Eq.  (26) . 


The  cost  of  repairing  a  crack,  C^(x) ,  depends  on  the  crack  size  x.  For  the 
crack  size  smaller  than  0.03",  it  is  only  necessary  to  ream  the  fastener  hole  to 
the  next  hole  size  and  hence  the  cost  of  repair,  C2(x),  is  identical  for  all 
x£0.03".  For  the  crack  size  greater  than  0.03",  a  regular  retrofit  procedure  may 
be  needed.  For  even  larger  cracks,  a  patching  replacement  of  the  entire  component 
may  be  necessary  which  is  much  more  expensive.  Therefore,  the  cost  of  repairing 
a  crack  size  x,  C^(x),  is  in  general  an  increasing  function  of  the  crack  size  x 
Likewise  C^(x)  may  also  be  a  discontinuous  function  of  x  as  described  above  [Ref.  10]. 


The  probability  of  detecting  a  crack  size  g.(x)dx  in  the  size  range  (x,x+dx) 
depends  on  the  resolution  capability  of  a  particular  NDI  technique  used  at  the  jth 
inspection  as  well  as  the  statistical  distribution  of  the  crack  size  at  T , .  It  is 
observed  from  Eqs.  < 22)  —  (.30)  that  the  estimation  of  the  crack  exceedance  and  the 
maintenance  cost  depends  on  p(i,t)  and  g,(x).  These  two  quantities,  along  with  the 
structural  reliability,  will  be  derived  ^in  the  next  section. 
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SECTION  VI 


DERIVATION  FOR  p(i,  t) 


AND 


g.(x) 


Even  under  well-controlled  laboratory  conditions,  the  detection  of  a  given 
crack  size  is  a  statistical  variable  and  it  should  be  specified  statistically. 

Let  F^(x)  be  the  probability  of  detecting  a  crack  size  x  and  F*(x)=l-F  (x)  be 
the  probability  of  missing  it.  Clearly,  both  F^(x)  and  F*(x)  depend  on  the 
crack  size  x  and  the  resolution  capability  of  a  particular  NDE  procudure  [e.g., 
Refs.  48,49].  Various  functional  forms  for  F^(x)  have  been  suggested  [11,16] . 

6.1  WITHOUT  INSPECTION  AND  REPAIR  MAINTENANCE 

Without  scheduled  inspection  and  repair  maintenance,  the  probability  p(i,T) 
that  the  crack  size  a(x)  of  a  detail  in  the  ith  stress  region  will  exceed  a  value 
x^  at  the  service  time  T  is  given  by 


p(i,T)  =  p[a(T)>x1]  =  1  -  Fa(T)(*1>  (31) 

Since  a(x)  is  related  to  the  initial  crack  size  a(0)  through  Eq.  (21)  in 
which  t  =0,  t^^x,  i.e., 

a(0)  *  W[a(x)  ,xl  (32) 


Eq.  (31)  becomes 

p(i,T)  =  p[a(0)>W(x1#T)]  =  P  [a (0)^(1)] 


or 


p(i#T)  =  1  -  Fa(0) Cy1(x)] 


(33) 


in  which  F 
given  by 


a  (0) 


(x)  is  the  distribution  function  of  the  equivalent  initial  flaw  size 
Eq.  (\2)  and  argument  y^(x)  is  related  to  x^  as 


y1(x)  «=  \Hxyt) 


(34) 
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The  value  of  y^(x)  given  by  Eq.  (34)  is  obtained  from  the  master  curve  for 
a  given  value  of  x^  as  shown  in  Fig.  (20). 

Thus  p(i,x)  is  obtained  by  substituting  Eq.  (12)  into  Eq.  (33)  as 
follows: 

a. 


p(i,T)  =  1  -  exp | 

=  0 


"L - J 


(35) 


in  which  x  is  given  by  Eq.  (13).  Note  that  p(i,x)  is  a  function  of  the  crack 
size  x^  when  Eq.  (34)  is  substituted  into  Eq.  (35). 

For  a  special  case  in  which  the  crack  growth  rate  in  the  ith  stress  region 
follows  Eq.  (8)  ,  i.e.,  da(t)/dt  =  Q.a\(t),  the  master  curve  can  be  written, 

after  integration,  as  follows 


a  (0)  = 


a  ( t) 


[1  +  aci<t)c^.t]1//ci 


=  W[a(t)  ,t] 


(36) 


where  c_.=b.-l.  Hence  y^(x)  appearing  in  Eq.  (35)  can  be  expressed  using  Eq.(36) 
in  the  following: 


yx(T)  = 


[l  +  xJcQj] 


I7c.  =  w[xi'x] 


(37) 


6.2  SCHEDULED  INSPECTION  AND  REPAIR  MAINTENANCE 


When  the  durability  critical  component  is  subjected  to  a  scheduled  inspection 
and  repair  maintenance  at  T^ ,  T  ,  T^,  ...with  service  intervals  x^ ,  x^,  ...as 

shown  in  Fig.  (1)  ,  the  solution  ' for  p(i,x)  is  more  involved.  In  the  first 

service  interval,  the  entire  population  of  the  initial  flaws  is  propagated  follow¬ 
ing  the  master  curve  as  shown  in  Fig.  (2  ).  During  the  first  inspection  and  repair 
performed  at  t=^>  some  of  the  cracked  details  will  be  detected  and  repaired.  The 
crack  size  of  the  repaired  details,  referred  to  as  the  renewal  details,  is  assumed 
to  have  the  same  distribution  as  that  of  the  EIFS  given  by  Eq.  (12).  in  other 
words,  the  possible  cold-work  effect  and  others  are  neglected.  Thus,  after  each 
inspection  and  repair  maintenance,  a  new  population  (renewal  details)  is  introduced 
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Figure  20 


Computation  of  y, (t)  and  y  Or) 
1  cr 
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with  certain  probability  depending  on  the  crack  size  distribution  prior  to  inspec¬ 
tion  as  well  as  the  detection  probability  F  (x)  of  a  particular  NDI  procedure. 
Furthermore,  the  statistical  distribution  of  the  crack  size  of  those  undetected 
details  is  modified  and  truncated  by  the  F*(x),  i.e.,  the  probability  of  nondetec¬ 
tion.  U 


Consequently,  the  distribution  of  the  crack  size  of  the  entire  population 
is  continuously  shifted  and  changed  (to  the  large  value)  in  service  due  to  fatigue 
crack  propagation.  The  distribution  is  also  subjected  to  truncation  and  modifi¬ 
cation  after  each  inspection  and  repair  maintenance.  Likewise,  since  an  additional 
new  population  is  introduced  after  each  inspection  and  repair  maintenance,  the 
statistical  distribution  of  the  crack  size  after  the  jth  inspection  and  repair 
maintenance  will  consist  of  j+1  different  populations  as  shown  in  Fig.  (21). 

Referring  to  Figure  (21),  f^^  ^ (x)  is  the  probability  density  of  all  cracks 

at  T  just  before  the  jth  inspection^and  repair  (see  Fig.  (1)).  F  (x)  is  the 
probability  of  being  able  to  detect  (and  repair)  cracks  of  size  x.  ^When  these 
two  are  multiplied,  the  product  represents  the  density  of  all  cracks  which  are 
found  and  repaired.  The  area,G(j),  under  the  curve  represents  the  percentage  of 
repaired  holes  that  will  have  a  probability  density  of  cracks  corresponding  to  that 
of  the  EIFS,  denoted  by  ^a(Q)(x)*  Cracks  which  are  missed  stay  at  the  same  length; 
that  is,  they  will  not  be  repaired.  This  portion  is  shown  as  fa(T) (x)  then  mul¬ 
tiplied  by  Fd(x)  where  Fg(x)=l-Fp(x)  is  the  probability  of  missing  a  crack  size  x. 
After  repair, the  probability  density  of  cracks,  fa(T.+)(x)>  i-s  composed  of  the  sum 

jf  ( jc )  F*  ( x)  3 

of  the  undetected  cracks,  a(T-i)  D  ’  and  repaired  cracks  f  .  .(x)G(j), 

a(0) 

The  inspection  and  repair  procedure  can  be  input  into  the  analysis  at  any  time 
desired.  Thus,  the  effect  of  different  inspection  or  repair  procedures  or  intervals 
can  be  ascertained. 

As  a  result,  a  new  population  is  introduced  after  each  inspection  and  repair 
maintenance.  The  variation  of  the  probability  density  of  the  flaw  size  in  service 
due  to  crack  propagation  and  maintenance  procedures  is  displayed  in  Fig.  (21). 

Let  ^a(o)^  ^e  the  probability  density  function  of  the  EIFS,  a(0).  Since 

the  probability  density  function  is  the  derivative  of  the  distribution  function, 
one  has 


'a  (0) 


(x)  = 


dF 


a  (0) 


(x) 


dx 


(38) 


in  which  the  distribution  function  F  , 

a  (0) 


(x) 


of  a(0)  is  given  by  Eq .  (12). 
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21  Effect  of  Repair  on  Crack  Population 


Substitution  of  Eq.  (12)  into  Eq.  (38)  leads  to  the  expression  for  the 
probability  density  function  of  the  equivalent  initial  flaw  size 


(39) 


The  probability  p(i,f)  that  the  crack  size  of  a  detail  in  the  ith  stress  region 
at  any  service  time  x=T,+t  will  exceed  a  value  x^  }  see  Fig.  (l)  t  for  j  =  0,1,2,... 
is  given  by  ^ 


P(i,T) 


fa(Tj+t)  ^  dx 


(40) 


in  which  f  +t)  (x)  t*ie  probability  density  function  of  the  crack  size  a(T.+t), 

and  t  is  any  value  smaller  than  or  equal  to  T.+^=  T.+^-T.  for  j=0,l,2,...,  i.e., 

t<T .  -T  .  Hence,  a(T  +t)  is  the  crack  size  at  any  service  time  T.+t  in  the  j+lth 
J+l  J  j  J 

service  interval  (see  Fig.  (1)). 

The  probability  density  function,  f  (x) ,  of  the  crack  size  a(T.+t),  at 

a(T^+t)  j 

any  service  time  T=T^.+t(j=0,l,2, .  .  .)  is  derived  in  the  Appendix  for  both  the  special 
case  and  the  general  use.  The  results  are  given  in  the  following; 


6.2.1  In  The  First  Service  Interval,  j=0 

In  the  first  service  interval,  p(i,i)  is  identical  to  that  given  by  Eqs.  (34) 
and  (35).  The  transformation  of  the  density  functions  with  the  relationship 
a(0)=W  [  a(t) , t  ]  can  be  found 


fa(t)(x)  *  =  fa(0) 


(41) 
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in  which 


yx(t) 

lx(t) 


W(x,t) 

dW(x,  t)  „ 
dx 


(42) 


where  the  crack  size,  y^t),  and  the  corresponding  crack  growth  rate,  dy^O/dt 

(slope),  are  obtained  from  the  crack  growth  master  curve  as  shown  in  Fig.  (22  (a)). 
The  crack  growth  rate,  dx/dt,  when  the  crack  size  is  x,  is  also  obtained  from  the 
master  curve.  Fig. (22  (a)). 


6.2.2  After  the  First  Service  Interval,  j=l,2. 


fa(T.+t)(x) 


—  A  .  , ,  (t;  j )  +  £  G(  j-k+1)  A.  (t;  j) 

D+1  k=l  * 


Vt;j>  “  |X  fa<0>  £yjt(tj3>  J1Jc<t'-3> 


(43) 


(44) 


k=2,3,..., j+1 


y  (t;j)  -  y1  t  +  £  xl  -  wlx,  t  +  £  xHj 

1L  n=  j-m+2  nJ  L  n=3-m+2  J 


dy  (t;  j) 

m _ _ 


dx 


Fdym(t;j)  ]  L 

y-^iA 


dx 

.dt. 


xn 


which  x  is  the  nth  service  interval  as  shown  in  Fig.  (1) . 
n 


(  (45) 


The  crack  size,  y  (t;j),  and  the  crack  growth  rates  (slopes),  dy^(t;j)/dt 
and  dx/dt,  appearing  in  Eqs.  (43)-<45)  are  determined  from  the  crack  growth 
master  curve  as  shown  in  Fig.  (22).  The  analytical  expressions  of  these  quanti¬ 
ties  for  a  special  case  are  given  in  the  Appendix.  Furthermore,  G(j)  (j=l,2,...) 
appearing  in  Eq.  (43)  represents  the  probability  of  detecting  a  crack  of  any  size 
at  the  jth  inspection. 


G(j) 


-  s: 


(X)  dx 


(46) 
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dx/d 


Figure  22  Utilization  of  Master  Curves 


5 


in  which  g^.(x)dx  is  the  probability  of  detecting  a  crack  in  the  size  range 


(x,x+dx)  at  the  jth  inspection  T 


9j(X>  "  fa(Tj)(x)FD(x) 


(47) 


Where  T  =0  and  f  .  (x)  =  f  ..  (x)  is  the  probability  density  function  of 

0  a(T.)  a(T  +TJ 

3  j-1  3 

the  crack  size  a(T.)  at  T.  prior  to  the  jth  inspection,  given  by  Eq.  (43). 


The  probability  of  detecting  (or  repairing)  a  crack  in  the  size  range 
(z  ,z  )  at  the  jth  inspection,  i.e.,  T.  (j=l,2, . . .)  flight  hours,  denoted  by 

q.(*l,«2).  follows  from  Eq.  (47)  as 


'j<Zl'22>  =  /  2 


(x)  dx 


(48) 


and  the  probability  of  detecting  a  crack  size  greater  than  z^  at  T.,  denoted  by 
qj  (Zf,00)  ,  is  therefore  ^ 


(x)  dx 


(49) 


which  is  referred  to  as  the  crack  exceedance  for  the  repaired  cracks  at  the  jth 
maintenance.  It  is  obvious  from  Eq.  (49)  .  that 


(x)  =  “dqj(Xjco)/dx 


(50) 
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SECTION  VII 


STRUCTURAL  RELIABILITY, 
CRACK  DETECTION  AND  REPAIR 


The  probability  of  detecting  a  crack  in  the  size  range  (x,x+dx)  at  the  jth 
inspection  is  given  by  Eq.  (47)  which  depends  on  the  detection  probability  F^(x) 
of  a  particular  NDI  technique.  Various  NDI  procedures  can  be  used  for  crack 
detection,  such  as  visual  inspection,  liquid  penetrant  inspection,  ultrasonic  inspec 
tion,  radiographic  inspection,  etc.  The  detection  probabilities  F^(x)  under  well- 

controlled  laboratory  conditions  as  a  function  of  the  crack  size  x  have  been  avail¬ 
able  for  various  NDI  techniques  [e.g..  Refs.  48-49].  The  experimental  results  have 
been  used  as  a  basis  to  characterize  the  functional  form  of  F^(x)  for  the  purpose 
of  engineering  analyses  [Refs.  11,14,16,17] .  One  functional  form  was  proposed 
[Refs.  11,14]. 


(51) 


in  which  a^  is  a  crack  size  below  which  the  crack  cannot  be  detected  and  a^  is  a 
crack  size  beyond  which  the  crack  can  be  detected  with  certainty.  Values-  a^,  a^  and 
m  should  be  determined  to  best  fit  the  experimental  results  of  a  particular  NDI  pro¬ 
cedure.  Equation  (51)  has  been  used  [e.g..  Refs.  11,13,14,15]  for  engineering 
analyses. 

Another  functional  form  for  the  crack  detection  probability  F^(x)  was  proposed 
in  Ref.  16  as  follows: 


Fd(x)  =  0 


for  x£a^ 


r  -fMx-a  n 
=  Cl[l  -  e  1  J  for  ax<x 


(52) 
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in  which  a^  has  the  same  meaning  as  that  in  Eq .  (51)  and  is  a  value  slightly 
less  than  unity  to  allow  for  a  small  probability  of  missing  a  big  crack  during 
inspection.  Again,  the  parameter  values  of  a  ,  c^  and  should  be  chosen  to  best 
fit  the  experimental  results. 

Sufficient  information  on  the  effectiveness  of  field  inspection  is  not  avail¬ 
able  for  a  characterization  of  F^Cx).  However,  Eqs.  (51)  and  (52)  can  be  used 
with  larger  values  of  a^,  a^,  and  c^  should  the  field  data  become  available.  It 
is  the  general  consensus  that  the  effectiveness  of  the  field  inspection  is  much 
lower  than  that  under  laboratory  conditions. 

If  a  denotes  the  critical  crack  size  associated  with  the  maximum  design  load 
cv 

P  then  the  failure  of  a  detail  (such  as  a  fastener  hole)  occurs  when  its  crack 
xx 

size  exceeds  a  .  Without  a  scheduled  inspection  and  repair  maintenance,  the  prob¬ 
er 

ability  of  failure  of  one  detail  in  the  ith  stress  region  of  the  component  within 
the  service  interval  (0,T),  denoted  by  P^O,!),  follows  from  Eq .  (31)  as 


P±(0,t)  =  P[a(x)>_a 


cr>  -  1  "  Pa(T)tacr> 


Following  the  similar  derivation  for  Eqs.  (33)  and  (34)  P^(0,t)  can  be 

shown  as 


P^O.T)  =  l  -  Fa(0)  [yor(T>] 


in  which  F  ^ (x)  is  the  distribution  function  of  the  equivalent  initial  flaw  size 
represented  by  Eq.  (12),  and  y  (T)  is  related  to  a^  through  Eq.  (21)  as 


ycr(T)  =  w(acr,t) 


For  a  given  value  of  the  critical  crack  size  a  ,  y  (t)  is  obtained  from 
°  cr  cr 

Eq.  (55)  using  the  crack  growth  master  curve  as  shown  in  Fig.  (20). 


The  probability  of  the  survival  of  all  the  N.  details  in  the  ith  stress  region 
of  the  component  in  the  service  interval  (0,t)  is 


r  iN- 

[l  -  Pi(0rT)J  x 
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Hence  the  probability  of  survival  of  the  entire  durability  critical  com¬ 
ponent  consisting  of  m  stress  regions  in  the  service  interval  (0,x)  denoted  by 
R(t),  is 


R(T) 


(57) 


and  -the  probability  of  failure  of  the  entire  component  in 


(0,t)  is 


Pf(T)  -  1  -  R(t)  =  1 


m  r 

n  |i-p. (o,t) 

i=l «-  1 


(58) 


Under  the  scheduled  inspection  and  repair  maintenance  as  shown  in  Fig.  (1) , 
the  probability  of  failure  of  a  detail  in  the  ith  stress  region  within  the  jth 
service  interval  (T.  ^,T.)  is  given  by 


pi‘Tj-i*V  = 


/ 


cr 


”a(Tj) 


(x)  dx 


for  1| 2/ • » • 


(59) 


in  which  f  (x)  =  f  ^  +T  ^ (x)  is  the  probability  density  function  of  the 

j  j — i  j 

crack  size  a(T^)  prior  to  the  jth  inspection  as  derived  in  Eq.  (43).  The  prob¬ 
ability  that  a  detail  in  the  ith  stress  region  will  survive  j  service  intervals 

(0,T.)  is 
J 


5 

J! 

1=1 


(60) 


in  which  TQ=0.  The  probability  that  the  ith  stress  region  consisting  of 
details,  will  survive  in  the  service  interval  (0,T  J  is 


f  3  r 

•l 

N. 

n  fi 

-  P.  (T0  . #T0 )  I 

JL 

1 1=1  L 

1  1-1'  1J 

(61) 
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Consequently,  the  probability  of  survival  of  the  entire  durability  critical 

component  which  consists  of  m  stress  regions  in  the  service  interval  (0,T  )  is 

3 


S(V  =  ilfilt1  •  Wi-v])"1 


(62) 


and  the  probability  of  failure  of  the  durability  critical  component  in  (0,T  )  is 

j 


■f'V  =1- j1(j1[1-pi!Vi'V]Ni) 


(63) 
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SECTION  VIII 


DEMONSTRATION  OF  METHODOLOGY 


In  order  to  demonstrate  the  applicability  of  the  durability  analysis 
methodology  developed  herein,  the  economic  life  of  a  hypothetical  durability 
critical  component  of  a  fighter  aircraft  is  considered  for  simplicity.  A 
scheduled  inspection  and  repair  maintenance  procedure  is  also  considered 
and  its  effect  on  the  economic  life  of  the  component  is  examined.  The 
component  consists  of  only  one  stress  region  having  100  identical  fastener 
holes  where  the  maximum  stress  level  of  the  fighter  spectrum  for  each 
fastener  hole  is  amax  =  ^  ksi.  Cracks  are  assumed  to  occur  at  the  fastener 
holes  only.  The  fastener  holes  of  the  component  is  assumed  to  be  identical 
to  those  of  the  WPF  data  set  presented  in  Table  1,  page  22. 


Test  data  of  TTCI  at  0.03  inches  for  WPF  and  XWPF  data  sets  were  pre¬ 
sented  in  Tables  1  and  2,  respectively,  and  summarized  in  Table  4.  The 
three-parameter  Weibull  distribution,  Eq.  (2),  was  used  to  best  fit  the 
experimental  data  as  shown  in  Figs.  (8)  and  (9).  The  parameter  values  of 
a,  8,  e  and  8y  for  y=95%  confidence  level  are  presented  in  Tables  3  and  5. 

The  crack  growth  master  curves  and  the  associated  crack  growth  rates  obtained 
using  the  CGR  computer  code  for  crack  growth  damage  accumulation  were  pre¬ 
sented  in  Figs.  (13)  to  (16).  The  corresponding  b  and  Q  values  in  the  small 
crack  size  range  are  given  in  Table  5. 


With  the  aid  of  Eq.  (9),  the  backward  extrapolations  have  been  carried 
out  to  obtain  the  corresponding  equivalent  initial  flaw  sizes  from  the  test 
data  of  times  to  crack  initiation.  The  results  are  presented  in  Table  4 
and  plotted  in  Figs.  (23)  and  (24)  as  circles.  Also  plotted  in  Figs.  (23) 
and  (24)  as  a  solid  curve  and  a  dashed  curve  are  the  theoretically-derived 
distributions  for  the  equivalent  initial  flaw  size  given  by  Eq.  (12).  The 
solid  and  dashed  curves  are  for  8  and  8y  for  y=  95%  confidence  level,  respec¬ 
tively. 


Considering  the  theoretical  EIFS  distribution  represented  by  the  solid 
curve  of  Fig.  (23)  and  the  crack  growth  master  curve  represented  by  curve  1 
of  Fig.  (13) ,  we  can  demonstrate  the  crack  damage  statistics  in  service  and 
the  economical  life  of  the  hypothetical  component  using  the  analytical 
methodology  developed  in  the  previous  sections „ 
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Table  4  TTCI  AND  EIFS  FOR  WPF  AND  XWPF  DATA  SETS 


WPF 

XWPF 

TTCI 

EIFS 

TTCI 

EIFS 

(Flight  Hours) 

(Mils) 

(Flight  Hours) 

(Mils) 

6560 

3.281 

6111 

0.903 

9310 

1.656 

6140 

0.893 

10631 

1.232 

7432 

0.566 

10888 

1.165 

7535 

0.547 

11348 

1.057 

7723 

0.514 

11630 

0.997 

8422 

0.412 

12051 

0.914 

8528 

0.399 

13373 

0.704 

8611 

0.389 

13540 

0.682 

8905 

0.357 

13753 

0.655 

8967 

0.350 

13765 

0.653 

9074 

0.339 

13783 

0.651 

9080 

0.339 

13936 

0.633 

9312 

0.317 

14106 

0.613 

9955 

0.265 

14124 

0.611 

10244 

0.245 

14129 

0.610 

10452 

0.231 

14256 

0.596 

10894 

0.206 

14347 

0.586 

11030 

0.199 

14400 

0.580 

11035 

0.199 

14433 

0.577 

11036 

0.199 

15492 

0.477 

11326 

0.185 

15600 

0.468 

11418 

0.181 

15807 

0.451 

11472 

0.178 

16079 

0.430 

11557 

0.175 

16087 

0.429 

11695 

0.169 

16719 

0.386 

11700 

0.169 

16741 

0.384 

11710 

0.168 

16752 

0.384 

11907 

0.160 

16979 

0.369 

12042 

0.155 

17043 

0.365 

12114 

0.153 

17056 

0.365 

12303 

0.146 

17396 

0.345 

12387 

0.143 

17556 

0.336 

12500 

0.139 

18139 

0.305 

12623 

0.136 

18696 

0.279 

12660 

0.134 

20140 

0.224 

13325 

0.116 

20449 

0.214 

16004 

0.067 

23767 

0.134 
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Without  inspection  and  repair  maintenance  procedures,  the  probability  p(i,x), 
that  a  detail  will  have  a  crack  size  larger  than  x  at  any  service  time  t  is  com¬ 
puted  from  Eq.  (35)  .  The  average  percentage  of  details  L(x)/N  ,  see  Eqs.  (26) 

and  (28)  ,  having  a  crack  size  larger  than  x^  is  plotted  in  Figure  (25)  as_a  A 

function  of  service  time  T  in  flight  hours.  In  Figure  (25)  >  the  ordinate  L(x)/N  , 
denoted  by  the  average  percentage  of  crack  exceedance,  is  the  average  number  of 
details,  while  the  abscissa  is  the  crack  size  x^.  For  instance,  at  8,000  flight 
hours,  the  average  percentage  of  fastener  holes  having  a  crack  size  exceeding  0.03 
and  0.05"  are  2%  and  0.8%,  respectively.  As  expected,  the  average  percentage  of 
crack  exceedance,  L(x)/N*,  increases  as  the  service  time  X  increases.  Curves  in 
Figure  (25)  are  referred  to  as  the  average  crack  exceedance  curves. 


Table  5  Model  Parameters  Used  for  Methodology 
Demonstration  for  EIFS 


WPF  * 

XWPF  ** 

b 

1.2165 

1.26 

Q 

%  • 

0.9247  x  10~3 

0.2328  x  10 "2 

a 

4.8634 

5.6995 

0 

14,957  hours 

11,372 

€ 

1,312 

0 

@7 

14 , 240 

10,864 

No  load  transfer 
15%  load  transfer 
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(X)  NOIlSNnd  NOlineiUiSiQ 


Figure  23  Cumulative  Probability  Versus  EIFS  for 
Weibull  Compatible  Distribution  (WPF) 
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DISTRIBUTION  FUNCTION  F 


1.0< 


Figure  24  Cumulative  Probability  Versus  EIFS  for 
Weibull  Compatible  Distribution  (XWPF) 


65 


As  discussed  in  Section  V,  the  number  of  details  having  a  crack  size  exceed¬ 
ing  at  any  service  time  T  is  a  statistical  variable  following  the  Binomial  dis¬ 
tribution  represented  by  Eq.  (22).  When  the  number  of  details  (fastener  holes) 
is  large,  the  Binomial  distribution  of  Eq.  (22)  can  be  approximated  by  the  normal 
distribution  with  the  mean  and  the  variance  given  by  Eqs.  (24)  and  (25),  respec¬ 
tively,  As  a  result,  Figure  (25),  represents  an  approximation  of  the  percentage 

of  fastener  holes  having  a  crack  size  exceeding  x  (an  abscissa)  with  50% 
probability.  1 

From  the  standpoint  of  the  durability  requirement  and  the  economic  life  cri¬ 
teria,  the  central  measure  of  the  crack  exceedance,  i.e.,  exceedance  with  50% 
probability,  as  shown  in  Figure  (25),  may  not  be  conservative  enough,  although 
it  provides  information  for  estimating  the  average  repair  cost  as  will  be  described 
later.  Hence,  it  is  desirable  to  compute  the  number  (or  percentage)  of  details 
having  a  crack  size  exceeding  x^  at  any  service  time  x  with  certain  probability, 
in  particular  the  low  probability  of  exceedance.  This  objective  can  be  accom¬ 
plished  using  either  the  Binomial  distribution  given  by  Eq .  (22)  or  the  normal 
distribution  for  approximation. 

For  simplicity  in  presenting  the  numerical  results,  the  service  time  of 
X=8000  flight  hours  is  considered.  The  percentage  of  cracks,  L(x)/N  ,  Eqs.  (26) 
and  (28)  ,  exceeding  the  crack  size  x^  (abscissa)  is  plotted  in  Figure  (26)  for 
various  exceedance  probabilities .  For  instance,  at  x=8,000  flight  hours,  the  prob¬ 
ability  is  0.1  that  3.764%  of  the  total  cracks  will  exceed  the  crack  size  0.03"  as 
indicated  by  an  open  circle.  Similarly,  the  probabilities  are  0.5  and  0.8,  respec¬ 
tively,  that  2%  and  0.803%  of  the  total  cracks  will  exceed  0.03"  (see  open  circles). 
It  is  further  observed  from  Figure  (26)  that  the  probabilities  are  0.1,  0.5  and 
0.8,  respectively,  that  1.97%,  0.8%  and  0.06%  of  the  total  cracks  will  exceed  0.05" 
as  shown  by  the  solid  circles.  Figure  (26)  indicates  that  as  the  probability  of 
exceedance  decreases,  the  percentage  of  cracks  exceeding  a' certain  crack  size  x^ 
increases  as  expected. 

The  numerical  results  presented  in  Figures  (25)  and  (26)  for  the  crack 
exceedance  with  certain  probability  does  not  account  for  the  confidence  level, 
since  the  value  of  8  representing  the  characteristic  time  to  crack  initiation  is 
used  in  the  analyses  (see  Table  5) .  When  a  confidence  level  Y  is  further  consid¬ 
ered,  then  $y  should  be  used.  The  value  of  8y  fox  y=95%  confidence  level  given 
in  Table  5  is  now  used  in  the  analyses.  Numerical  results  are  presented  in 
Figures  (27)  and  (28)  which  correspond  to  a  95%  confidence  level.  For  instance, 
it  is  observed  from  Figure  (28)  that  at  8,000  flight  hours,  the  probability  is 
0.1  that  4.51%  of  the  total  cracks  will  exceed  0.03"  with  95%  confidence  level  as 
indicated  by  a  circle.  In  other  words,  at  8,000  flight  hours,  4.51%  of  the  total 
cracks  will  exceed  0.03"  with  0.1  probability  and  95%  confidence. 
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AVERAGE  PERCENTAGE  OF 
CRACK  EXCEEDANCE 


Figure  25  Crack  Exceedance  Curves  for  50%  Confidence  Level 


AT  8000  HOURS 


CRACK  SIZE  (MILS) 


26  Crack  Exceedance  and  Reliability 
for  507o  Confidence  Level 


AVERAGE  PERCENTAGE  OF 


(957.  CONFIDENCE) 


CRACK  SIZE  (MILS) 


Figure  27  Crack  Exceedance  Curves  (95%  Confidence) 
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A  comparison  between  Figs.  (25)  and  (27)  indicates  that  the  average  per¬ 
centage  of  cracks  exceeding  any  given  crack  size  x^  increases  when  a  higher  con¬ 
fidence  level  is  considered.  Further  comparison  between  Figs.  (26)  and  (28) 
shows  that  the  percentage  of  cracks  exceeding  any  given  crack  size  x^  with  any 
given  probability  increases  as  a  higher  level  of  confidence  is  taken  into  arrnunt 

The  economic  repair  crack  size,  ae,  is  arbitrarily  defined  as  the  max¬ 
imum  crack  size  below  which  the  least  expensive  repair  procedure  can  be  used 
i.e.,  reaming  the  fastener  hole  to  the  next  hole  size.  For  the  above  defini¬ 
tion,  ae  is  usually  between  0.03"  and  0.05"  depending  on  the  location  and  the 
fastener  hole  size.  If  0.03"  is  used  as  the  economic  repair  crack  size,  then 
the  average  percentage  of  details  having  a  crack  exceeding  ae=0.03"  is  plotted 
as  a  function  of  the  service  time  as  shown  by  Curve  1  of  Figure  (29).  This 
curve,  in  fact,  can  be  constructed  from  Figure  (25)  by  drawing  a  vertical  line 
passing  through  0.03".  The  average  percentage  of  cracks  exceeding  ae  is  also 
the  percentage  of  cracks  exceeding  a£  with  0.5  probability.  Hence,  Curve  1  is 
referred  to  as  the  exceedance  curve  for  the  economic  repair  crack  size  a£  with 
probability  0.5.  .The  exceedance  curves  for  the  economic  repair  crack  size  with 
probability  5%  and  10%,  respectively,  are  also  computed  and  plotted  in  Figure 
(29)  as  dashed  curves.  It  is  observed  from  the  figure  that  Curve  1  increases 
rapidly  after  8,000  flight  hours.  If  5%  of  the  cracked  details  exceeding  :0.03" 
with  probability  50%  is:  considered  as  an  economical  limit,  then  the  possible 
economic  life  is  observed  to  be  9,400  flight  hours. 

When  the  economic  repair  crack  size  a  is  0.05",  one  can  construct  in  a  similar 
fashion  the  exceedance  curve  with  probabilities  50%  (Curve  1) ,  10%  (dashed  curve) 
and  5%  (dashed  curve)  as  shown  in  Figure  (30).  The  possible  economic  life  in  this 
case  is  found  to  be  10,500  flight  hours  if  5%  of  the  cracks  exceeding  0.05"  with 
probability  50%  is  considered  as  an  economic  limit. 

In  the  same  manner,  one  can  also  construct  the  exceedance  curves  for  the 
economic  repair  crack  size  a  =0.03"  with  probabilities  0.5  (solid  curve),  0.1 

P 

(dashed  curve)  and  0.05  (dashed  curve)  as  a  function  of  service  life  when  a  con¬ 
fidence  level  of  y=95%  is  taken  into  account.  The  results  are  displayed  in  Figure 
(31).. 


In  order  to  demonstrate  the  effect  of  inspection  and  repair  maintenance  pro¬ 
cedures  on  the  economic  life,  an  NDI  procedure  having  a  high  resolution  capability 
is  assumed.  The  functional  form  for  the  detection  probability  F^(x)  given  by  Eq . 

(51)  is  used  with  a^=0.01",  a^^O.l"  and  m=0.5.  Such  a  detection  probability  F^(x) 
is  plotted  in  Figure  (32)  as  a  function  of  the  crack  size  x  . 

The  cost  formulation  for  the  economic  life  criteria  as  represented  by  Eq.  (30) 
is  based  on  the  average  (expected)  maintenance  cost,  including  the  average  costs 
of  inspection  and  repair.  Since  the  average  cost  of  repair  is  related  to  the  average 
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PERCENTAGE  OF  CRACKS 
EXCEEDING  aQ 


5  6  8  10  12  14  16 

FLIGHT  HOURS,  103 


Figure  29  Crack  Exceedance  Versus  Service  Time  Showing  the 
Effects  of  Repairs  (ag  =  0.03") 
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PERCENT  OE  CRACKS 


Figure  30 


Crack  Exceedance  Versus  Service  Time  Showing  the  Effects 
of  Repairs  (ag  ■  0.05") 
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Figure  31  Crack  Exceedance  Versus  Service  Time  (95%  Confidence) 
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number  of  detected  cracks  in  various  size  ranges,  g.(x)dx,  the  effect  of  the 
maintenance  procedures  on  the  economic  life  will  be^ investigated  on  the  basis 
of  the  average  crack  exceedance  curve  (i.e.,  with  50%  exceedance  probability) 
as  shown  by  Curve  1  of  Figures  (29)  and  (30)  with  50%  confidence  level. 

For  the  case  where  a  =0.03",  Figure  (29)  ,  the  first  inspection  and  repair 

e 

is  performed  at  T^=9500  flight  hours.  The  resulting  average  exceedance  curve 

(with  50%  probability)  is  obtained  by  use  of  Eqs.  (43)  to  (45)  and  plotted 
in  Figure  (29)  as  Curve  2.  Owing  to  the  repair  procedures,  there  is  a  sudden 
drop  from  (a)  to  (^)  (i.e.,  Curve  1  to  Curve  2)  at  as  indicated  in  Figure 

(29)  •  According  to  Curve  2  of  Figure  (29)  ,  the  possible  economic  life  (5% 

allowable  exceedance  over  a  with  50%  probability  as  defined  previously)  is  now 

6 

extended  to  11,000  flight  hours. 


If  the  second  inspection  and  repair  is  further  performed  at  T  =11,000  flight 
hours,  then  the  resulting  average  exceedance  curve  is  computed  using  Eqs.  (43) 
to  (45)  and  plotted  in  Figure  (29)  as  Curve  3.  The  economic  life  is  thus 
extended  to  12,000  flight  hours.  In  a  similar  manner,  additional  inspections  and 
repair  maintenance  procedures  are  performed,  respectively,  at  T^=12,000  and  T^= 
13,000  flight  hours.  The  resulting  average  crack  exceedance  curves  are  plotted 
in  Figure  (29)  as  Curves  4  and  5,  respectively.  It  is  observed  from  Figure  (29) 
that  there  is  a  sudden  drop  of  the  average  percentage  of  crack  exceedance  at  each 
inspection  and  repair  maintenance,  which  is  exactly  the  average  percentage  of 
cracks  repaired  (or  detected)  during  each  maintenance. 


Thus,  if  the  durability  critical  component  is  subjected  to  a  scheduled  inspec¬ 
tion  and  repair  maintenance  procedure  at  T^,  T^ ,  T^,  and  T^,  then  the  average  crack 
exceedance  curve  will  follow  Curve  1  up  to  ©  and  then  follows  the  path  ©-*€>* 

Curve  5 . 

For  the  case  where  the  economic  repair  crack  size  is  0.05",  the  inspection 

and  repair  maintenance  procedures  are  performed  at  T^=10,500,  1^=12,500,  T^= 

14,000  and  T  =15,250  flight  hours,  respectively.  The  resulting  average  crack 
4 

exceedance  curves  are  presented  in  Figure  (30).  If,  for  instance,  the  inspec¬ 
tion  and  repair  maintenance  are  performed  at  T^  and  T^  only,  then  the  resulting 
average  crack  exceedance  curve  will  follow  Curve  1  up  to  and  then  follows 

the  path  ©— ©— ®-©  — ►Curve  3. 


The  average  percentage  of  details  having  a  crack  size  larger  than  z^  which 
are  repaired  at  T . ,  denoted  by  q  .  (z^,<»)  (see  Eq .  (49)),  is  plotted  in  Figure  (33) 

for  j=l,2,3,  and  4  for  the  case  where  a  =0.03".  Curves  in  Figure  (33)  are 
referred  to  as  the  average  exceedance  curves  for  the  repaired  cracks.  For  instance, 
one  percent  of  the  details  having  a  crack  size  larger  than  0.1"  is  expected  to  be 
repaired  at  T^=9,500  flight  hours  as  indicated  by  a  circle  on  Curve  1  of  Figure 
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AVERAGE  PERCENTAGE  OT 
REPAIRED  CRACK  EXCEEDANCE 


Figure  33  Average  Exceedance  Curves  for  Repaired  Details 


(33)  •  The  exceedance  curves  depicted  in  Figure  (33)  indicate,  on  the  aver¬ 

age,  the  number  of  details  as  well  as  their  corresponding  crack  size  to  be 
repaired  during  each  inspection  and  repair  maintenance.  These  curves  contain 
the  information  needed  for  the  computation  of  the  average  repair  cost  repre¬ 
sented  by  Eq.  (30)  ,  since  g .  (x)=-dq  .  (x,°°)  /dx,  see  Eq  ,  (50). 

It  is  observed  from  Figures  (29)  and  (30)  that  the  effect  of  inspection 

and  repair  maintenance  procedures  on  the  economic  life  is  small  even  using  a 

high  resolution  NDI  procedure.  It  is  found  that  for  other  NDI  techniques  with 

poor  detection  capability,  the  effect  is  practically  minimal.  The  reasons  are 

given  as  follows:  (i)  To  extend  the  economic  life,  the  NDI  procedure  should  be 

able  to  detect  with  high  probability  the  crack  size  below  the  economic  repair 

crack  size  0.03"  in  order  to  prevent  cracks  from  propagating  beyond  a  .  Current 

6 

NDI  techniques  may  not  achieve  such  a  capability  with  a  comparable  cost,  and  (ii) 
As  soon  as  the  crack  exceedance  curve  over  a  starts  to  increase  rapidly  [see 
Figures  (29)  and  (30)  the  majority  of  the  crack  population  is  already  in  the 
vicinity  of  a  ,  and  hence  a  significant  extension  of  the  economic  life  implies 

G 

an  extensive  repair  such  that  the  cost  of  repair  may  be  uneconomical. 
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SECTION 


I  X 

CONCLUSIONS  AND  DISCUSSIONS 


A  durability  analysis  methodology  for  aircraft  structures  has  been  devel¬ 
oped.  The  durability  analysis  is  based  on  the  statistical  approach  which  allows 
for  the  determination  of  a  possible  economic  life  using  one  of  the  following 
criteria:  (i)  the  percentage  of  crack  population  exceeding  the  economic  repair 
crack  size  with  certain  probability  and  confident  levels,  (ii)  the  ratio  of  the 
cost  of  maintenance,  including  the  costs  of  inspection  and  repair,  to  the  initial 
cost.  Although  the  inspection  and  repair  maintenance  procedures  have  a  signifi¬ 
cant  impact  on  the  reliability  and  safety  of  aircraft  structures,  its  effect  in 
extending  the  economic  life  is  shown  to  be  limited. 

In  the  present  durability  analysis,  the  input  initial  fatigue  quality  can 
either  be  the  time  to  crack  initiation  (TTCI)  or  the  equivalent  initial  flaw 
size  (EIFS) .  The  three-parameter  Weibull  distribution  is  used  for  the  statisti¬ 
cal  distribution  of  the  TTCI.  The  statistical  distribution,  which  is  derived 
from  that  of  the  TTCI  using  a  particular  crack  propagation  law,  is  suggested 
for  the  distribution  of  the  EIFS.  Such  a  derived  distribution  possesses  the 
physical  meaning  for  the  fatigue  wear-out  process  and  it  is  shown  to  fit  the 
EIFS  better  than  other  distributions. 

The  CGR  computer  code  is  used  to  calculate  the  master  curve  (crack  growth 
damage)  in  each  stress  region  of  the  entire  durability  critical  component  in  ser¬ 
vice.  The  binomial  distribution  is  used  to  compute  the  statistical  distribution 
of  the  number  of  fastener  holes  having  a  crack  in  any  size  range.  The  crack 
exceedance  with  certain  probability  and  confidence  is  then  estimated  as  a  func¬ 
tion  of  the  service  time.  The  relaibility  of  the  durability  critical  component 
in  service  is  also  obtained.  A  scheduled  (nonperiodic)  inspection  and  repair 
maintenance  procedure  is  considered,  in  which  the  statistical  uncertainty  of  the 
NDI  procedure  and  crack  growth  damage  of  the  repaired  details  are  taken  into 
account.  The  average  cost  of  maintenance,  including  the  costs  of  inspection  and 
repair,  is  then  formulated  as  a  function  of  service  time. 

While  the  second  criterion  for  determining  the  possible  economic  life  has 
been  formulated  and  relevant  statistics,  such  as  the  crack  exceedance  for  the 
repaired  cracks,  have  been  obtained,  it  is  not  demonstrated  numerically.  This  is 
due  to  a  lack  of  realistic  information  of  the  field  inspection  cost  as  well  as 
the  cost  of  repairing  one  crack  in  a  certain  size  range.  It  appears  that  the 
first  criterion  of  the  economic  life  may  be  more  appropriate  at  the  preliminary 
design  stage  for  the  durability  requirements. 
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Although  the  present  approach  and  formulation  are  demonstrated  using  the 
example  of  fastener  holes,  they  are  applicable  to  cracks  which  occur  at  radii 
or  shear  webs  with  holes,  as  long  as  the  crack  growth  master  curve  is  given  or 
computed . 

There  are  circumstances  where  the  usage  or  missions  of  aircraft  may  deviate 

significantly  at  certain  service  time  of  its  design  life.  Such  a  usage  change 

can  easily  be  taken  care  of  in  the  present  durability  analysis.  For  instance, 

when  the  usage  change  occurs  at  t  flight  hours,  then  the  new  crack  growth  master 

curve,  which  corresponds  to  the  new  usage,  should  be  used  in  the  analysis  after 

t  and  all  the  analysis  procedures  remain  the  same, 
c 

Although  the  reliability  of  the  durability  critical  component  has  been  formu¬ 
lated  within  the  framework  of  the  durability  analysis,  it  is  not  demonstrated 
numerically.  It  will  be  investigated  in  the  Phase  II  effort  whether  the  relia¬ 
bility  problem  is  important  in  the  durability  requirement. 

In  the  present  approach,  the  crack  growth  master  curve  under  spectrum  loading 
is  considered  deterministic,  i.e.,  the  statistical  variability  of  crack  growth 
parameters  is  neglected.  This  may  be  a  reasonable  approximation,  since  the  statis¬ 
tical  variability  of  the  crack  growth  parameters  is  smaller  than  that  of  the  equiva¬ 
lent  initial  flaw  size,  and  hence  the  statistical  variability  of  the  crack  exceed¬ 
ance  curves  (or  economic  life)  is  essentially  due  to  that  of  the  EIFS.  Such  an 
approximation  simplifies  the  entire  analysis  tremendously.  The  effect  of  the  sta¬ 
tistical  variability  of  the  crack  growth  damage  on  the  economic  life  of  aircraft 
structures  will  be  investigated  in  the  Phase  II  effort. 

It  should  be  emphasized  that  the  statistical  variability  of  the  crack  growth 
damage  is  important  and  should  be  accounted  for  in  the  damage  tolerant  analysis. 

In  the  damage  tolerant  analysis,  the  initial  flaw  size  is  specified  and,  hence, 
the  statistical  uncertainty  in  structural  life  prediction  is  essentially  con¬ 
tributed  to  by  the  uncertainty  of  the  crack  growth  damage  when  the  service  loads 
are  idealized  as  spectrum  loadings. 
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APPENDIX 


DERIVATION  OF  PROBABILITY 

DENSITY  FUNCTION  f  .  (x) 

a(T.+t) 


(A)  SPECIAL  CASE 

In  order  to  facilitate  the  derivation,  the  crack  growth  rate  equation  repre¬ 
sented  by  Eq.  (8)  will  first  be  used  for  simplicity.  Eq.  (8)  will  only  serve 
as  a  convenient  vehicle  to  explain  the  mathematical  derivation.  The  general 
results  without  the  limitation  of  Eq.  (8)  will  be  described  later. 

Integration  of  Eq.  (8)  from  t^  to  t^  yields 


a(t2) 


a(t^  °  it;  a°(t2)c0(t2-t1)]vc  ■  (A-i, 


or 

a(t,) 

a(t2)  =  1 - c - ^7T7H  =  w 

[1  -  a^t^cQ^-t^]1'0  11  2  1J 

in  which  t^St^,  W  ^[.]  is  the  inversed  function  of  W  [. ]  and 


(A— 2} 


C  =  b  -  1 


(A- 3) 


For  t^  =0  and  t^  =t,  Eqs .  (A-l)  and  (A-2)  reduce  to 


a  (0)  = 


a  (t) 


[1  +  aC(t)cQt]1/c 


=  W[a(t) ,t] 


a(t)  =  7 - iTT - 7T7Z  "  W^IatOJ^t] 

[1  -  ac(0)CQt]1/c 


(A-4) 
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where  a(0)  is  the  equivalent  initial  flaw  size. 

Eqs .  (A-l)  to  (A-4)  will  serve  as  the  transfer  functions  to  transform  the  dis¬ 
tribution  function  of  the  EIFS  into  that  of  the  flaw  size  at  any  service  time. 

It  should  be  mentioned  that  the  values  of  Q  and  b  (or  c=b-l)  vary  from 
one  stress  region  to  another.  However,  the  values  of  Q  and  b  (or  c=b-l)  appear¬ 
ing  in  the  distribution  of  EIFS  are  obtained  from  one  data  set  of  coupon  specimens 
under  the  same  loading  spectra  but  at  a  particular  maximum  stress  level.  As  a 
result,  Q  and  c  appearing  in  this  appendix  should  be  written  as  Q.  and  c.  except 
Eqs.  (A-7)  and  (A-8) .  Because  of  simplicity  in  notation,  the  subscript  i  has  been 

dropped  and  they  should  not  be  confused  with  Q  and  c  for  the  EIFS. 

(1)  First  Service  Interval  (0,T^) 

In  the  first  service  interval,  the  entire  initial  crack  population  is  sub¬ 
jected  to  propagation  following  Eq .  (A-4). 

Let  F  .  .  (x)  and  f  .  .  (x)  be  the  cumulative  distribution  function  and  the 
a (t)  a(t) 

probability  density  function,  respectively,  of  the  crack  size  a(t)  at  t  flight 

hours  in  the  first  service  interval.  Then  both  F  ,  .  (x)  and  f  ,  .  (x)  can  be 

a(t)  a (t) 

derived  from  those  of  the  initial  flaw  size  F  (x)  and  f  (x) ,  respectively, 

a(0)  a(0) 

through  the  transformation  of  Eq .  (A-4)  as  follows: 


pa(t)(x)  =  PU(t)<x]  =  Pta(0)<yi(t)]  =  Fa(0)  [yi(tn 

or 

Fa(t)  (x)  =  Fa(0)  [yl(t)  J  (A-5) 


in  which  y^(t)  is  a  value  of  a(0)  corresponding  to  a  value  x  of  a(t).  They  are 
related  through  Eq .  (A-4)  by  replacing  a(0)  and  a(t),  respectively,  by  y^(t)  and 
x. 


yx(t) 


_ X 

[1  +  XCcQt]1/C 


W(x,t) 


Substituting  Eq. 


(12)  into  Eq .  (A-5),  one  obtains 


(A- 6) 
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Fa(t) 


(x) 


-  1 


I 


YlUIlXu  (A-7) 


in  which  y  (t)  is  a  function  of  x  given  by  Eq .  (A-6)  and  x  is  given  by  Eq.  (13). 
1  u 

The  probability  that  the  crack  size,  in  the  ith  stress  region  at  T =t  ser¬ 
vice  hours,  will  exceed  x^  is  given  by 


Pu,T)  =  i  -  ra(T)  (Xl> 


=  1  -  exp 


[V^C  ( t )  -aQC-cQe*| 

p  L  coe  J  ; 


yi(T>ixu 


(A- 8) 


=  0 


t  yx^)>\ 


in  which  Eq.  (A-7)  has  been  used  and  y  (x)  is  related  to  x  through  Eq .  (A-4) 
as  11 


yx(t)  = 


[l  +  x^cQx] 


T7c 


=  Wlx. ,t] 


(A-9) 


As  expected,  p(i,t)  obtained  in  Eq.  (A-8)  is  identical  to  that  derived  in 
the  text  in  which  no  inspection  and  repair  maintenance  is  scheduled. 

The  transformation  of  the  probability  density  function  from  f^^(x)  to 
fa(t)(x)  can  be  found  in  any  text  book  on  probability  and  statistics. 


aCt) 


(x)  =  f. 


(0)  [y1(t)]j1<« 


(A— 10) 


in  which  y^(t)  =  W(x,t)  is  related  to  x  through  Eq.  (A-6)  and 


J1(t)  = 


dy^t) 

dx 


_  dW(x,t) 
dx 


_l+xccQt_ 


where  Eq.  (A-6)  has  been  used. 


(A-ll) 
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At  the  end  of  the  first  service  interval,  i.e.,  t  =  =  T^,  the  probability 

density  function  of  the  crack  size  a(T^)  given  by  Eq .  (A-10)  becomes 


;a(T,)(x)  =  fa(0) 


(A-12) 


in  which  it  follows  from  Eqs.  (A-6)  and  (A-ll)  that 


yi,Ti>  = 


Ji(Ti>  = 


[1+xccQt1] 

1 


I7c 


=  W{x,T,) 


r  c  i*+1 

|_1  +  x  cQiJ  • 


dW(x,T1) 

dx 


(A-13) 


The  component  is  then  subjected  to  inspections.  The  probability,  g^(x)dx, 
of  detecting  a  crack  in  the  size  range  (x,x+dx)  at  t=  =T^  is  therefore 

gi(x)dx  =  FD(x)fa(  j (x)dx  (A-14) 

in  which  f  ,  (x)  is  given  by  Eq.  (A-12)  and  F_(x)  is  the  probability  of  detect- 

aCT^)  D 

ing  a  crack  size  x. 

The  probability  of  detecting  a  crack  in  the  size  range  z and  z^  during  the 
first  inspection,  denoted  by  q  (z  ^2) 

f* 2  f*2 

ql(zl,z2)  =  J  g^xjdx  =  J  FD(x)fa{T  j  (x)  dx  (A-15) 
21  Z1  1 

and  the  probability  of  detecting  a  crack  of  any  size  at  T^,  indicated  by  G(l),  is 


/CO  /•  oo 

(x)  dx  =  J  FD(x)  fa(T^j  (x)  dx  (A-16) 
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(2)  Second  Service  Interval  (T^.T^) 


It  is  assumed  that  when  a  crack  is  detected,  it  is  repaired.  Thus,  after 
the  first  inspection  and  repair  performed  at  T  ,  the  details  (fastener  holes) 
consist  of  two  different  populations;  (1)  the  population  where  cracks  have  been 
detected  and  repaired,  referred  to  as  the  renewal  details,  and  (ii)  the  popula¬ 
tion  where  the  cracks  have  not  been  detected  and  hence  details  have  not  been 
repaired.  It  is  further  assumed  that  the  crack  size  distribution  of  the  renewal 
details  is  identical  to  that  of  the  equivalent  initial  flaw  size  a(0) .  Conse¬ 
quently,  the  probability  density  function  of  the  crack  size  a(T+)  right  after 

the  first  inspection  and  repair,  denoted  by  f  ,  .  (x) ,  consists  of  two  parts 

a(TJ  ) 


fa(T*)  (X) 


G(1) fa(0) *x) 


+  FD(x)fa(Tl)  (X) 


(A-17) 


in  which 


F*(x)  =  1  -  pD(x) 


(A-18) 


is  the  probability  of  missing  a  crack  size  x  during  inspection,  and  G(l)  and 
fa(T  are  given  by  Eqs*  (A~16)  and  (A-12) ,  respectively. 

It  is  obvious  that  the  first  term  in  Eq.  (A-17)  represents  the  contribution 
from  the  renewal  (repaired)  population  with  probability  G(l) ,  while  the  second 
term  is  contributed  by  the  old  population  (which  survives  the  inspection  with 
probability  F*(x)). 


Substitution  of  Eq.  (A-13)  into  Eq .  (A-12)  and  then  into  Eq .  (A-17)  yields 


the  following  expression  for  f  (x) 

a(Tp 


W)<*>  -  c«1»fa(0)(x>  +  FS(*>fa(o)r - 2 - 175] 

L(1+xccQt1)  J 


[i+xccQt  J 


+1 


(A- 19) 

In  the  second  service  interval  (T  ,T  ) ,  the  two  different  crack  populations 
given  by  Eq.  (A-19)  are  subjected  to  propagation  in  service.  Let  the  probability 
density  function  of  the  crack  size  a(x)  at  any  service  time  x  =T  +t  in  (T^jT^)  be 
denoted  by  (x)  .  Then,  f^^  +(_^  (x)  can  be  obtained  from  fa(T+)  (x) 
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given  by  Eq .  (A-19)  through  the  transformation  of  Eq .  (A-l)  with  t^T-j+t, 
t1=^1  [Ti!=t1»  see  Fig-  CD  ]  »  i-e-j 


a(T,) 


afT^+t) 


|^l+ac(T1+t)  cQtJ 


X7c  =  W  [a(i^+t) ,  t] 


(A-20) 


Following  the  same  procedures  described  above  by:  (i)  replacing  the  argument 
of  x  in  Eq.  (A-19)  by  y  (t)  =  W(x,t)  ^Eq .  (A-6)]  and  (ii)  multiplying  Eq .  (A-19) 
by  J^(t),  one  obtains  the  following  expression  for  f^^  (x)  , 


^(Tj+t)  (x>  Gl(1)  fa(0)  Iyl(t) 

+  Fgly^t)  ]fa(o)  ]J2(i1+t)J1(t) 


(A-21) 


in  which  it  can  be  shown  that 


y^T^+t)  =  w(x,Tl+t)  = 


[l+XCCQ(T]L  +  t)J 


1/c 


J2(TJL+t)  =  =  ri±XCcQt 


r 

~dx  I 


I  +1 


1+X  CQft+Tj^Jj 


(A-22) 


At  the  end  of  the  second  service  interval,  i.e.,  T  =T  +T2  (°r  t=T2)>  the 
probability  density  function  of  the  crack  size  a^^)  f°ll°ws  from  Eq .  (A-21)  as 


fa(T0)<X)  “  G'1>fa(0)‘5,l(T2,1Jl(T2» 

1  2  (A-23) 

+  lfa(0)  >J2l  WW 


The  probability  of  detecting  a  crack  in  the  size  range  (x,x+dx)  at  T 
denoted  by  g^MdXjis 


2’ 
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g2(x)dx  =  FD(x)fa(T^)  (x) 


(A-24) 


in  which  f  (x)  is  given  by  Eq .  (A-23) .  The  probability  of  detecting  a 
i 2' 


crack  size  between  z^  and  z^  is 


q2(z 


l'z2>  "  J 


92(x) 


dx 


rD 


tx)fa(T2)(x)dx 


(A-25) 


and  the  probability  of  detecting  a  crack  of  any  size  during  the  second  inspec¬ 
tion  performed  at  is 


G(2)  =  q2(0,»)  = 


/CD 

FD(x)fa(T2) 


(x)  dx 


(A-26) 


(3)  Third  Service  Interval  (T  ,T^) 

The  probability  density  function  of  the  crack  size  a(T+)  right  after  the 
second  inspection  and  repair  consists  of  an  additional  contribution  from  the 
renewal  (repaired)  population  with  probability  G(2)  (i.e.,  produced  by  the 
second  inspection  and  repair) , 


fa(T2)  *x)  =  Gt2)fa(0)*x)  +  FD(x)  fa(T2)  (x)  (A-27) 

in  which  f  .  (x)  is  given  by  Eq.  (A-23). 
a(i2; 

Substitution  of  Eq.  (A-23)  into  Eq.  (A-27)  yields 

fa(Th(x)  -  G(2)fa(0)(x)  +  G(1)F£(x>fa(0)  ‘yi(T2),J1(T2) 

(A-28) 

+  F*(x)  F*  [y^(i2)  ^  J2  ^Tl+T2^  Jl^x2^ 

The  crack  size  a(T2+t)  at  any  service  time  '1r2+t(t£  To)  in  the  third  service 
interval  is  related  to  the  crack  size  a(T+)  through  Eq .  (A-l)  in  which  t^=T2 , 
T2+t ,  i .  e . , 

a{T  +t) 

a(T2)  =  7— - TI7H  =  Wta(T2+t)  ,t]  (A-29) 

I 1+a  (T2+t)cQt| 
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The  probability  density  function  of  a(T  +t) ,  denoted  by  f  .  (x)  is 

2  a (T  +t) 


obtained  from  f 


a(T+) 


(x)  through  the  transformation  of  Eq .  (A-29)  employing  the 


same  procedures  described  above  by:  (i)  substituting  Eqs.  (A-22)  with  t=  into 
Eq.  (A-28) ,  (ii)  replacing  x  in  Eq.  (A-28)  by  y  (t)=W(x,t)  given  by  Eq .  (A-6) 
and  (iii)  multiplying  Eq .  (A-28)  by  J  (t);  with  the  results 


a(T2+t) 


(x)  =  G ( 2)  f 


a  (0) 


ty^ ( t)  3 J-^t) 


G(l)F5[yi(t)]fa(0) 


[yi(r2+t) ]J2(T2+t) J1(t)  +  F* ty1(t) ]F*[yi(t2+t) ] 

fa(0)  [yl(Tl+T2+t)  JJ2(Ti+T2+t)  J2(x2+t)  Jl(t)  <A“3°) 

in  which  it  can  be  shown  that 


J2*Tl+T2+t*  = 


1+x  cQ(i2+t) 

L l+xCcQ(x1+T2+t) 


l+1 


(A-31) 


and  yi(t1+T2+t)  is  8iven  by  Eq .  (A-22)  in  which  the  argument  T^+t  is  replaced 

by  T1+T2+t' 

At  the  end  of  the  third  service  interval  T^T^+T^,  *'e,>  t=T3>  ^q .  (A-30) 
becomes 


fa(T3)  “  G*23  fa(0)  *Y1*T33  3  J1*T33  +  G<  ^  F£  *  T33  3  fa  (0) 

ty1(T2+T3)  ]J2(t2+T3)  J1(t3)  +  F*[y]L(T3)  ]F*[yi(T2+T33 

fa (0)  fy1(x1+x2+T3)3J2(Tl'fT2+T3)J2(T2+T3)Jl(T3) 

(A-32) 


The  probability  of  detecting  a  crack  in  the  size  range  (x,x+dx)  at  the 
third  inspection  is 


Mix  =  FD(x>fa(T3) 


(x)dx 


(A-33) 


and  the  probability  of  detecting  a  crack  of  any  size  at  the  third  inspection  is 


G  ( 3) 


g3(x)dx  = 


/ 


FDix)fa(T3)  <x)dx 


(A-34) 
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in  which  f  /rn  .  (x)  is  given  by  Eq.  (A-32)  . 

a(V 


(4)  ith  Service  Interval  (T.,T..,) 

- j -  j  j+i 

In  a  similar  fashion,  the  general  expressions  for  the  probability  density 
function  of  the  crack  size  a(T.+t)  in  the  jth  service  interval  where  tiT_.+^-T^ 
(j>l)  can  be  shown  as  "* 

*  ,  j 

fa(T.+t)(x)  “  Aj+l(t; j)  +  z  G(j+l-k)Ak(t;j)  for  j=l,2,... 

3  k“1  (A-35) 


in  which 


* 


k-1 

"  FDtyn.(t;j,) 
vm— l  J 

fa(0) tyk(t;j>3 

(A— 36) 


y2(t;j)  =  yx(t)  = 


x 


(l+xccQt) 


"I7c  =  W(x'fc> 


y_(t;  3) 
m 


=  y ,  ft  +  i  T  1  =  — - - - . - 

lL  nJ  i  x  if 

l  L  n=j~m+2  nJj 


(A— 37) 


=  W  x,  t  +  Z  T  1 

L  n=j-m+2  nJ 


Jx(t;j)  = 


Jx(t)  = 


dyx(t) 

dx 


[l+xccQt] 


i+1 


J2(t; j) 


l+xccQt 


H  +1 
i  c 


1+X  CQ(t+Tj) 


(A-38) 


J  (tt  j) 
m 


’ 

l+xCcQ 

r*  +  i  o 

n= j-m+3  _ 

j 

l+xccQ 

t  +  l  Tn 

L  n= j-m+2  -J 

J 

l+1 


m=3. 4f 


i 


*  Term  Equal  to  1.0  for  k  =  1 
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It  can  be  easily  shown  from  Eq.  (A-38)  that 


=  ",  = 
m=i 


l+xccQ 


f  3 

t  +  Z 
L  n=j- 


= j-k+2 


IjL  ( t ;  j )  =  J^tyj) 


Consequently,  Eq .  (A-36)  can  be  simplified  as 
A1  (t;j)  =  fa(o)  ^yi  (t;j)  ^  11  (tiJ) 

(k-1  ) 


-  -L  ± 

A^(t.-j)  =  Pn  Fg[ym(t;j)]|  fa(0)  [yk<t;jl]  ik(t;j) 

m=l  ' 


(A-39) 


(A-40) 


in  which  I,  (t;j)  is  given  by  Eq .  (A-39)  and  y  (t,j)  is  given  by  Eq.  (A-37) . 
k  m 

(B)  GENERAL  CASE 

For  the  general  situation,  the  crack  growth  rate  represented  by  Eq .  (8) 
may  not  be  applicable  because  of  the  crack  geometry  and  other  factors,  especially 
in  the  large  crack  size  regime.  As  a  result,  a  general  crack  growth  computer  pro¬ 
gram,  such  as  CGR,  should  be  used  to  generate  the  crack  growth  master  curve.  Then, 
the  crack  growth  damage  accumulation  can  be  expressed  by  Eqs .  (4-1)  or  (4-2)  (also 
see  Eqs.  (A-l)  and  (A- 2) ,  i.e., 

aft^}  =W[a(t2)# 

or 

a(t2)  =  w"1[a(t1)/  t2-tx] 

in  which  t^-t^. 

For  t^=0  and  t2=t,  Eqs.  (A-41)  and  (A-42)  becomes 

a (0)  =  W  [a (t)  ,  t]  (A-43) 


(A-41) 


(A-42) 


a (t)  =  W**1  (a(0)  ,  t] 


(A-44) 


90 


Now,  Eqs.  (A-41)  to  (A-44)  will  serve  as  the  transfer  functions  to  transform 
the  statistical  distribution  of  the  initial  flaw  size  a(0)  into  that  of  the  flaw 
size  at  any  service  time. 

(1)  In  the  First  Service  Interval  (0,T^) 


Following  the  same  derivations  presented  in  Eqs.  (A-5)  to 
show  that  p(i,x)  is  given  by  Eq.  (35)  or  Eqs.  (A-8)  and  (A-9) 


(A-9) ,  one  can  easily 
in  which 


y1(r)  =  WfXj^T) 


(A-45) 


The  probability  density  function  f  (x)  of  the  crack  size  at  t  flight  hours 
is  obtained  from  °f  EIFS  through  the  transformation  of  Eq.  (A-43) (see 
Eq.  (A-10»  as 


fa(t)  (x)  =  fa(0)  Iyl(t) 


(A-46) 


in  which  f  ... 

a(0) 


(x)  is  given  by  Eq.  (39)  and 


yx(t)  =  W(x,t) 

-  [m]/  [if] 


(A-47) 


Note  that  Eqs.  (A-46)  and  (A-47)  are  identical  to  Eqs.  (41)  and  (42)  ,  respec¬ 

tively. 


(2) 


In  the  j+1  th  Service  Interval  (T. 


Using  the  transfer  functions  given  by  Eqs.  (A-41)  to  (A-44)  instead  of  Eqs. 
(A-l)  to  (A-4)  (for  the  special  case)  and  following  exactly  the  same  procedure 
derived  above  for  the  special  case,  one  can  show  that  the  probability  density  func¬ 
tion  f  .  .  (x)  of  the  crack  size  a(T.+t)  at  the  service  time  T,+t  in  the  j+1  th 

3  WjT'tJ  j  j 

service  interval  is  given  by  Eqs.  (43)  to  (47)  in  the  text. 
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The  solutions  given  by  Eqs.  (43)  to  (47)  can  also  be  obtained  without  going 
through  the  derivation  described  above  using  the  following  rationale. 

Referring  to  Fig.  (21),  the  probability  density  function  at  T.+t  flight  hours 
in  the  j+1  th  service  interval  consists  of  j  different  populations.  This  is  because 
an  additional  new  population  is  introduced  after  each  inspection  and  repair.  The 
first  tern  of  Eq .  (43)  , 


Vi(t)  = 


D 

n 

m=l 


fa(0)tyj+l't!3,)ID+l(t;j>  (A-48) 


is  the  contribution  from  the  original  population  starting  at  the  zero  service  time. 

The  term  f  .  ,[y  (t; j  )]  I  (t;j)  is  the  probability  density  function  of  the  flaw 

a (.U;1-  j+1  j+i 


size  a (T.+t)  directly  obtained  from  that  of  a(0)  through  the  transformation  of  Eq . 

(A-43)  in  which  t  is  replaced  by  T.+t.  However,  the  original  population  is  subjected 
to  j  times  of  inspections  and  repairs.  Only  the  fraction  of  this  population  which  is 
not  detected  during  each  inspection  has  a  contribution  to  the  density  function  of 
a(T_j+t).  The  probability  (or  fraction)  of  not  being  detected  during  j  times  of  inspec¬ 


tions  is  given  by 


(  j 

( m=l 


,  in  which  a  crack  size  x  at  the  m  th  inspection  cor¬ 


responds  to  a  crack  size  y^(t;j)  prior  to  service  (see  Eq .  (43)). 


The  last  term  of  Eq.(43)for  k= j ,  i.e.,  G(l)A,(t;j)  represents  the  contribution 

from  the  new  population  introduced  after  the  first  repair  performed  at  T^  (see  Fig.(l)) 
with  probability  G(l).  Similarly,  the  second  term  of  Eq .  (43)  for  k=l 


GC j) Ax(ty  j) 


(A-49) 


represents  the  contribution  from  the  new  population  introduced  after  the  j  th  repair 
performed  at  T..  The  fraction  of  cracked  details  repaired  at  T.  is  G(j).  In  Eq . 

(A-49),  A^(t;j)  is  given  by  Eqs.  (41)  and  (42)  which  indicates  the  density  func¬ 
tion  of  the  crack  size  at  T.+t  when  its  crack  size  at  T.  is  a(0)  (cracks  repaired  at 


J  J 

Tj).  is  obtained  from  the  probability  density  function  of  a(0)  through  the 

transformation  of  Eq .  (A-41)  with  t ,=T,  and  t  =T  +t,  i.e., 

1  j  2  3 


a (Tj)  =  W [a (Tj+t) ,  tj 


(A-50) 


92 


REFERENCES 


1.  MIL-STD-1530A,  "Aircraft  Structural  Integrity  Program,"  Dec.  1975. 

2.  MIL-A-8866B,  "Airplane  Strength,  Rigidity  and  Reliability  Requirements;  Repeated 
Loads  and  Fatigue,"  Aug.  1975. 

3.  MIL-A-8867B,  "Airplane  Structural  Ground  Tests,"  Jan.  1975. 

4.  C.  F,  Tiffany,  et  al.,  "Analysis  of  USAF  Aircraft  Structural  Durability  and 
Damage  Tolerance,"  AIAA  Structural  Durability  and  Damage  Tolerance  Workshop, 
Bethesda,  April  6-7,  1978. 

5.  M.  D.  Coffin  and  C.  R.  Tiffany,  "New  Air  Force  Requirements  for  Structural  Safety, 
Durability,  and  Life  Management,"  Journal  of  Aircraft,  Vol.  13,  No.  2,  1976, 

pp.  93-98. 

6.  H.  A.  Wood,  R.  M.  Engle,  J.  Gallagher  and  J.  M.  Potter,  "Current  Practice  on 
Estimating  Crack  Growth  Damage  Accumulation  with  Specific  Application  to  Struc¬ 
tural  Safety,  Durability  and  Reliability,"  Air  Force  Flight  Dynamic  Lab., 
AFFDL-TR-75-32 ,  WPAFB ,  1976. 

7.  H.  A.  Wood,  "The  Use  of  Fracture  Mechanics  Principles  in  the  Design  and  Analysis 
of  Damage  Tolerant  Aircraft  Structures,"  Fatigue  Life  Prediction  for  Aircraft 
Structures  and  Materials,  AGARD-LS-62,  1973,  pp .  4. 1-4. 3. 

8.  H.  A.  Wood,  "Structural  Integrity  Technology  for  Aerospace  Applications," 
Structural  Integrity  Technology,  edited  by  J.  P.  Gallagher  and  T.  W.  Crooker, 

ASME,  May  1979,  pp.  23-28. 

9.  S.  D.  Manning,  J.  N.  Yang,  M.  Shinozuka,  W.  R.  Garver,  et  al.,  "Durability 
Methods  Development,  Volume  I  -  Phase  I  Summary,"  Air  Force  Flight  Dynamics  Lab., 
AFFDL-TR-7 9-3118 ,  Sept.  1979. 

10.  J.  N.  Yang,  "Statistical  Estimation  of  Service  Cracks  and  Maintenance  Cost 

for  Aircraft  Structures,"  Journal  of  Aircraft,  Vol.  13,  No.  12,  1976,  pp .  929-937. 

11.  J.  N.  Yang  and  W.  J.  Trapp,  "Reliability  Analysis  of  Aircraft  Structures  Under 
Random  Loading  and  Periodic  Inspection,"  AIAA  Journal,  Vol.  12,  No.  12,  1974, 
pp.  1623-1630. 

12.  J.  N.  Yang  and  W.  J.  Trapp,  "Inspection  Frequency  Optimization  for  Aircraft 

Structures  Based  on  Reliability  Analysis,"  Journal  of  Aircraft,  Vol.  12,  No.  5, 
1975,  pp.  494-496.  ~  ’  . 

13.  J.  N.  Yang,  "Statistical  Approach  to  Fatigue  and  Fracture  Including  Maintenance 
Procedures,"  in  Fracture  Mechanics,  Univ.  of  Virginia  Press,  Sept.  1978,  pp .  559- 
577,  Proc.  of  2nd  International  Conf.  on  Fracture  Mechanics,  Washington,  D.C. 


93 


14.  J.  N.  Yang  and  W.  J.  Trapp,  "Reliability  Analysis  of  Fatigue-Sensitive  Air¬ 
craft  Structures  Under  Random  Loading  and  Periodic  Inspections,"  Air  Force 
Materials  Laboratory,  AFML-TR-75-29,  WPAFB,  1975. 

15.  M.  Shinozuka,  "Development  of  Reliability-based  Aircraft  Safety  Criteria: 

An  Impact  Analysis,"  AFFDL-TR-76-29 ,  WPAFB,  1977. 

16.  J.  R.  Davidson,  "Reliability  After  Inspection,"  Fatigue  of  Composite  Materials, 
ASTM-STP  569,  1975,  pp .  323-334. 

17.  J.  R.  Davidson,  "Reliability  and  Structural  Integrity,"  NASA  TM  X71934,  1973. 

18.  J.  L.  Rudd  and  T.  D.  Gray,  "Equivalent  Initial  Quality  Method,"  Air  Force 
Flight  Dynamic  Lab.,  AFFDL-TM-76-83 ,  WPAFB,  1976. 

19.  J.  L.  Rudd  and  T.  D.  Gray,  "Quantification  of  Fastener  Hole  Quality,"  Proc. 

18th  AIAA/ASME/SAE  Structures,  Structural  Dynamics  and  Materials  Conf . ,  1977  . 

20.  J.  N.  Yang,  "Statistical  Estimation  of  Economic  Life  for  Aircraft  Structures," 
Proc.  AIAA/ASME/ASCE/ASH  20th  Structures,  Structural  Dynamics  and  Materials  Conf., 
St.  Louis,  April  1979,  AIAA  paper  No.  79-0761,  pp .  240-248;  Journal  of  Aircraft, 
AIAA,  Vol .  17,  No.  7,  1980,  pp .  528-535. 

21.  J.  N.  Yang  and  S.  D.  Manning,  "Statistical  Distribution  of  Equivalent  Initial 
Flaw  Size,"  1980  Proceedings  Annual  Reliability  and  Maintainability  Symposium, 

San  Francisco,  Jan.  1980,  pp.  112-120. 

22.  S.  Eggwertz,  "Investigation  of  Fatigue  Life  and  Residual  Strength  of  Wing  Panel 
for  Reliability  Purposes,"  Probabilistic  Aspects  of  Fatigue,  ASTM-STP  511,  1972, 
pp.  75-105. 

23.  P.  R.  Abelkis,  "Fatigue  Strength  Design  and  Analysis  of  Aircraft  Structures, 

Part  I:  Scatter  Factor  and  Design  Charts,"  U.S.  Air  Force  Flight  Dynamic  Lab., 
AFFDL-TR-66-197 ,  WPAFB,  1967. 

24.  A.  M.  Freudenthal,  "Reliability  Analysis  Based  on  Time  to  First  Failure," 

Fifth  ICAF  Symposium,  Australia,  1967. 

25.  A.  M.  Freudenthal  and  G.  I.  Schueller,  "Scatter  Factors  and  Reliability  of  Air¬ 
craft  Structures,"  NASA  Grant  NGR-09-010-058 ,  The  George  Washington  University, 
Washington,  D.C. 

26.  A.  M.  Freudenthal,  "Reliability  Assessment  of  Aircraft  Structures  Based  on 
Probabilistic  Interpretation  of  the  Scatter  Factor,"  U.S.  Air  Force  Materials 
Lab.,  AFML-TR-74-198,  WPAFB,  1975. 


94 


27.  I.  C.  Whittaker  and  P.  M.  Besuner,  "A  Reliability  Analysis  Approach  to  Fatigue 
Life  Variability  of  Aircraft  Structures,"  U.S.  Air  Force  Materials  Lab., 
AFML-TR-69-65,  WPAFB,  1969. 

28.  I.  C.  Whittaker,  "Development  of  Titanium  and  Steel  Fatigue  Variability  Model 
for  Application  on  Reliability  Analysis  Approach  to  Aircraft  Structures,"  U.S. 

Air  Force  Materials  Lab.,  AFML-TR-72-236 ,  WPAFB,  1972. 

29.  J.  P.  Butler  and  D.  A.  Rees,  "Development  of  Statistical  Fatigue  Failure  Char¬ 
acteristics  of  0.125  in.  2024-T3  Aluminum  Under  Simulated  Flight-by-Flight 
Loading,"  U.S.  Air  Force  Materials  Laboratory,  AFML-TR-124,  WPAFB,  1974. 

30.  C.  S.  Sarphie  and  R.  S.  Watson,  "Evaluation  of  Reliability  Analysis  Approach  to 
Fatigue  Life  Variability  of  Aircraft  Structures  Using  C-130  In-Service  Operational 
Data,"  AFML-TR-70-272 ,  WPAFB,  1971. 

31.  C.  S.  Sarphie,  "Examination  of  Values  of  Reliability  Distribution  Parameters 
Using  Structural  Fatigue  Data  from  C-130  Aircraft  Service  Utilization,"  AFML- 
TR-70-272,  WPAFB,  Vol .  II,  1973. 

32.  L.  F.  Impellizzeri ,  et  al . ,  "Evaluation  of  Structural  Reliability  Analysis 
Procedures  as  Applied  to  a  Fighter  Aircraft,"  U.S.  Air  Force  Materials  Lab., 
AFML-TR-73-150 ,  WPAFB,  1973. 

33.  W.  S,  Johnson,  R.  A.  Heller,  and  J.  N.  Yang,  "Flight  Inspection  Data,  Crack 
Initiation  Times,  and  Initial  Crack  Size,"  Proc.  1977  Annual  Reliability  and 
Maintainability  Symposium,  Jan.  1977. 

34.  J.  N.  Yang  and  W.  J.  Trapp,  "Joint  Aircraft  Loading/Structures  Response  Statistics 
of  Time  to  Service  Crack  Initiation,"  Journal  of  Aircraft,  AIAA,  Vol.  13,  No.  4, 
1976,  pp.  270-278. 

35.  P.  J.  Noronha,  et  al.,  "Fastener  Hole  Quality,  Vol.  I  &  II,  Air  Force  Flight 
Dynamics  Lab.,  AFFDL-TR-78-209 ,  WPAFB,  1978. 

36.  M.  Shinozuka,  "Durability  Methods  Development,  Volume  IV  -  Initial  Quality 
Representation,"  Air  Force  Flight  Dynamics  Lab.,  AFFDL-TR-79-3118,  Sept.  1979. 

37.  J.  P.  Gallagher  and  H.  D.  Stalnaker, "Predicting  Flight-by-Flight  Fatigue  Crack 
Growth  Rates,"  Journal  of  Aircraft,  Vol.  12,  1975,  pp .  699-705. 

38.  J.  P.  Gallagher  and  H.  D.  Stalnaker, "Methods  of  Analyzing  Fatigue  Crack  Growth 
Rate  Behavior  Associated  with  Flight-by-Flight  Loading,"  AIAA  paper  No.  74-367, 
presented  at  AIAA/ASME/SAE  15th  SDM  Conf.,  1974. 

39.  J.  P.  Gallagher,  "Estimating  Fatigue  Crack  Lives  for  Aircraft:  Techniques," 
Experimental  Mech.,  Vol.  16,  No.  11,  1976,  pp .  425-433. 


95 


40.  J.  P.  Gallagher,  "A  Generalized  Development  of  Yield  Zone  Models,"  Air  Force 
Flight  Dynamic  Lab.,  AFFDL-TM-74-27,  WPAFB,  1974. 

41.  J.  P.  Gallagher  and  T.  F.  Hughes,  "Influence  of  Yield  Strength  on  Overload 
Affected  Fatigue  Crack  Growth  Behavior  in  4340  Steel,"  AFFDL-TR-74-28 ,  WPAFB, 

1974. 

42.  G.  J.  Petrak  and  J.  P.  Gallagher,  "Predictions  of  the  Effect  of  Yield  Strength 
on  Fatigue  Crack  Growth  Retardation  in  Hp-9Ni-4Co-30C  Steel,"  ASME  paper  no. 
75-Mat-N . 

43.  R.  J.  Engle,  Jr.,  "Cracks  A  Fortran  IV  Digital  Computer  Program  for  Crack  Propa¬ 
gation  Analysis,"  AFFDL-TR-70-107 ,  WPAFB,  1970. 

44.  R.  M.  Engle  and  J.  L.  Rudd,  "Analysis  of  Crack  Propagation  Under  Variable 
Amplitude  Loading  Using  the  Willengborg  Retardation  Model,"  AIAA  paper  No. 

74-369,  15th  AIAA/ ASME/ AE  SDMConf.,  1974. 

45.  T.  R.  Porter,  "Method  of  Analysis  and  Prediction  of  Variable  Amplitude  Fatigue 
Crack  Growth,"  Engineering  Fracture  Mechanics,  Vol .  4,  1972. 

46.  0.  E.  Wheeler,  "Spectrum  Loading  and  Crack  Growth,"  ASTM  paper  No.  71-Met-X,  1971. 

47.  W.  Elber,  "Equivalent  Constant-Amplitude  Concept  for  Crack  Growth  Under  Spectrum 
Loading,"  ASTM  STP  595,  1976,  pp .  236-250. 

48.  P.  F.  Packman,  et  al.,  "Definition  of  Fatigue  Cracks  Through  Non-destructive 
Testing,"  Journal  of  Materials,  Vol.  4,  No.  3,  1969,  pp .  666-700. 

49.  B.  G.  W.  Yee,  et  al.,  "Assessment  of  NDE  Reliability  Data,"  NASA  contractor 
report,  NASA,  CR-134991,  1976. 


96 


